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We review methods developed in the gauge-string duality to treat energy loss by en- 
ergetic probes of a strongly coupled thermal medium. After introducing the black hole 
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density and the Poynting vector of the dual gauge theory. We summarize some efforts 
to compare these calculations to heavy ion phenomenology. 
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1. Introduction 

A long-standing hope, as yet unrealized, is that quantum chromodynamics (QCD) 



will be reformulated as a string theory. The gauge-string dualitj-^^^^ provides the 
closest approach to that goal so far attained. It provides useful computational meth- 
ods for studying strongly coupled gauge theories. The theory that is most accessible 
via these methods is Af — A super- Yang- Mills theory (SYM) in the limit of a large 
number of colors and large 't Hooft coupling Q Aspects of the progress in using the 
gauge-string duality to understand QCD have recently been reviewed at a pedes- 
trian level.* The aim of the current article is a more focused review of efforts to 

^ A/" = 4 super- Yang-Mills theory is a gauge theory whose matter content consists of gluons, four 
Majorana fermions in the adjoint representation of the gauge group, and six real scalars, also in 
the adjoint representation. All the fields are related to one another by the A/" = 4 supersymmetry, 
which completely fixes the Lagrangian once the gauge group and gauge coupling are chosen. Our 
interest is in the gauge group SU{N). 
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understand energy loss by energetic probes of a thermal medium in a strongly cou- 
pled gauge theory, such as SYM, which has a string theory dual; and to review how 
energy loss in SYM can be compared to energy loss in QCD. 

At least in simple cases like SYM, the response of an infinite, static, strongly 
coupled thermal medium to an energetic probe can be presumed to be hydrody- 
namical far from the probe, because hydrodynamic perturbations are the only long- 
wavelength modes available. In an infinite, interacting thermal medium, there is no 
radiation, because there are no asymptotic states. In a strongly coupled medium, 
it is not clear that there is a gauge-invariant distinction between coUisional energy 
loss and radiative energy loss. So the main questions are: 

(1) What is the rate of energy loss from an energetic probe? 

(2) What is the hydrodynamical response far from the energetic probe? 

(3) What gauge-invariant information can be extracted using the gauge-string du- 
ality about the non-hydrodynamic region near the probe? 

(4) Do the rate and pattern of energy loss have some meaningful connection to 
heavy ion phenomenology? 

In section [2] we briefly review the dual description of the thermal state of SYM 
as an AdSs-Schwarzschild black hole. The reader interested in a more extensive dis- 
cussion of the AdS/CFT duality is referred to various reviews in the literatureEEQ 
In section [3] we explain how to describe heavy quarks in A/" = 4 SYM using strings 
in AdS5, and in section [4] we extract the drag force acting on the quark via an 
"obvious" and "alternative" identification of parameters between SYM and QCD. 
Also in section 4, we consider how the string theory estimates of drag force relate 
to the measured nuclear modification factor for heavy quarks. 

The response of the stress-tensor to the motion of the quark is dual to the metric 
perturbations around an AdS black hole. These are studied in section |5] where we 
also discuss how the metric perturbations map into the stress tensor of the plasma. 
In section [6] we provide analytic approximations to the stress tensor both near 
to the moving quark and far from it. The full numerical solution is described in 
section [7] and its application to heavy-ion phenomenology can be found in section [8] 
In [Appendix A| we provide a glossary of mathematical notations used in the main 
text. 

2. The thermal medium as a black hole 

J\f — 4 SU{N) super- Yang-Mills theory at finite temperature can be described in 
terms of N D3-branes near extremality.'^ D3-branes are 3+1 dimensional objects 
on which strings may endE!i5l Each string can end on one of the N D3-branes, 
which eventually gives rise to the SU{N) gauge symmetry of SYM.^^ Having a 
non-vanishing tension, D-branes themselves warp the spacetime around them. The 
near horizon geometry produced by the D3-branes is AdSs-Schwarzschild times a 
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five-sphere, threaded by N units of flux of the self-dual Ramond-Ramond five-form 
of type IIB supergravity|i2l For most calculations of interest to us in this review, 
one can ignore the five-sphere and work just with the five non-compact dimensions. 

Our starting point is the Einstein equations of type IIB supergravity in the 
non-compact directions, which can be recovered from the action 



12 



(1) 



Here k is the five-dimensional gravitational coupling, G — detG^i^, and L is the 
radius of the S^. Standard relations based ultimately on the quantized charge of 
the D3-brane lead to 



L3 / n 
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.2 
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AdSs-Schwarzschild is a solution of the Einstein equations following from Q. Its 
line element is 



ds^ = Gn^dx'^dx'' = -h(z)dt^ + -I- -r-- , (3) 

where the "blackening function" is given by 

h{z) = 1 ~ 4 • (4) 

In this coordinate system, the conformal boundary of AdSs is located at z = 0. 
The AdSs-Schwarzschild solution has a horizon at z = z^f, whose temperature 

is 

. (5) 

According to the gauge-string duality, the temperature of the horizon ([5| is also 
the temperature of the thermal medium in the dual gauge theory!^ This medium is 
infinite and static. Its energy density e equals the mass per unit coordinate volume, 
d'^x, of the black hole, and its pressure p can also be straightforwardly computed 
as minus the free energy of the black hole. To leading order in the number of colors 
N and the 't Hooft coupling A, one finds 

These relations reflect the conformal invariance oi M — ^ super- Yang-Mills theory, 
which is exact even at flnite N and A. If we send z// — > oo, or equivalently T — ^ 0, 
we end up with a pure AdSs geometry and no black hole. 

Perturbations of the AdS5-Schwarzschild black hole with wavelengths much 
longer than l/T are described by relativistic fluid dynamics 1133 and the viscosity 
is known to be remarkably small; 

s 47r' ^> 
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again to leading order in the limit of large N and A. 

In summary: A thermal medium of SU{N) M — A super- Yang- Mills theory 
can be represented as a black hole in AdSs in the limit of large N and large 't 
Hooft coupling. Relations between the radius of AdSs and the five-dimensional 
gravitational coupling, between the temperature and the position of the horizon, 
between the energy density and the temperature, and between the viscosity and the 
entropy density, can all be derived starting from ten-dimensional type IIB string 
theory. 

3. The trailing string 

M — A super- Yang-Mills theory has no fundamentally charged quarks. Instead, 
its field content is the gluon and its superpartners under A/" = 4 supersymmetry: 
namely, four Majorana fermions and six real scalars, all in the adjoint representation 
of SU{N). The AdSs description of SYM makes no direct reference to these colored 
dynamical fields. The magic of AdS/CFT is to replace the strong coupling dynamics 
of open strings, whose low-energy quanta participate in the gauge theory, with the 
gravitational dynamics of closed strings (gravitons, for example) in a weakly curved 
background. For example, as we saw in section [2] a finite temperature bath in the 
gauge theory is replaced by a black hole in AdSs. The absence of colored degrees 
of freedom in the gravitational theory makes it less than obvious how to discuss 
energetic colored probes of the medium. 

A clue comes from the treatment of Wilson loops in the gauge-string duality 
A static, infinitely massive, fundamentally charged quark can be represented as a 
string hanging straight down from the boundary of AdSs, at z = 0, into the horizon. 
There is no contradiction with the previous statement that J\f — 4 super- Yang- Mills 
has no fundamentally charged quarks, because the quark is an external probe of the 
theory, not part of the theory. An anti-fundamentally charged quark is represented 
in the same way as a quark, except that the string runs the other way. (In type IIB 
string theory, strings are oriented, so this statement makes sense.) If a quark and an 
anti-quark are both present, it is possible for the strings running down into AdSs to 
lower their total energy by connecting into a shape like a catenary. From the point 
of view of the boundary theory, the string configuration gives rise to an attractive 
potential with a 1/r dependence, as conformal invariance says it must. 

Let's go back to a single isolated quark — which makes sense in A/" = 4 even at 
zero temperature because there is no confinement. The string dual to the quark can 
be described by its embedding in AdSs. The two-dimensional spacetime manifold 
swept out by the string is called a worldsheet and can be parameterized by two 
coordinates. If we use t and z to parameterize the string worldsheet, then the string 
in AdSs, or ™ AdSs-Schwarzschild, dual to a static quark at x — 0, is described by 
the equation 

static quark: x{t, z) — . (8) 
This is a solution to the string equations of motion with boundary conditions such 
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that the string endpoint on the boundary is stationary. By symmetry, ([8]) is the 
only possible solution with these boundary conditions. In the absence of a medium, 
it must be that the description of a moving quark can be obtained by performing a 
Lorentz boost on the description of a static quark. If the boost is in the xi direction, 
and we continue to use t and z to parameterize the worldsheet, then the equations 
describing the string dual to a quark moving in the absence of a medium are 

quark moving in vacuum: xi{t, z) — vt , xj_{t, z) = , (9) 

where v is the velocity of the quark and x± — (a;2, 0:3 )j^ Given that ^ solves 
the equations of motion of classical string theory, it's guaranteed that ^ does too, 
because the equations are invariant under the boost in the xi direction, and the line 
element ^ with /i = 1 is invariant under any boost acting on the (i, x) coordinates. 

When h 1, signaling the presence of a medium, a moving quark is not the 
same as a static one. Naively, one might nevertheless try to represent the moving 
quark in terms of the string shape described in Let's see why this is problematic. 
String dynamics is defined in terms of the metric on the worldsheet. We chose to 
parameterize the worldsheet by t and z, but any coordinates a" — ((t^,(T^) could 
have been chosen. For a general embedding x'' = x^{a^, cr^) of a string worldsheet 
into spacetime, the worldsheet metric is 

dx^dx^ 



What ( 10 1 says is that times and distances along the string are measured the 
same way as in the ambient spacetime. Using cr" = {t,z) and the ansatz ([9]), one 
immediately finds 

ff"/3 = TT n 1 • (11) 



This shows that the z direction on the worldsheet is always spacelike and that the 
t direction is timelike only when h > v^. Using Q, this condition is equivalent to 



z < z, = zhV^-v^ . (12) 

The part of the string worldsheet with z > is purely spacelike, which means 
that the string is locally moving faster than the speed of light. This does not make 
sense if we are aiming to describe the classical dynamics of a string moving in real 
Minkowski time. 



There is a simpler way to arrive at ( 12 1: any trajectory of a point particle in the 



bulk of AdSs-Schwarzschild with xi — vt must be spacelike — that is, the speed of 
the particle will exceed the local speed of light — if the inequality z < is violated. 



''It would be more proper to use {x^,x^,x^) in place of {xi, X2, X3), and reserve the notation 
Xfj, for Gfiux" . However, it simplifies notation to set Xi = 3;' for i = 1,2,3, and we will do this 
consistently. 
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This conclusion holds even when there are other components of the velocity, but 
the inequality is sharp only in the case where there aren't. 



Physically, what we learn from (12 1 is that it is harder and harder to move 
forward as one approaches the horizon. This evokes the idea that there must be 
some drag force from the medium. But how does one get at that drag force? An 
answer was provided by two groups jilli^ and closely related work on fluctuations 
appeared at the same time.li^ The string does not hang straight down from the 
quark: rather, it trails out behind it. If the shape is assumed not to change as 
the quark moves forward, and if it respects the SO(2) symmetry rotating the x± 
coordinates, then it must be specified by a small variant of ([9]): 

xi{t,z) =vt + ^{z) (13) 

for some function ^{z). If we insist that the quark's location on the boundary is 
xi — vt, then we must have ^ ^ as z — )■ 0. 

To determine ^{z), one must resort to the classical equations of motion for the 
string. These follow from the action 

1 



-'string 



d'cT^^, (14) 



27ra' 

where g — det g^js . The parameter a' is related to L and the 't Hooft coupling by 

= . (15) 
Here we define A = gYyiN, so that a quantity analogous to the coupling Us in QCD 



is q;ym = X/4:TtN. The classical equations of motion following from the action (13) 
take the form 

V„p"^ = , (16) 

where 

P%^-^g''^G^,d0x'' (17) 
is the momentum current on the worldsheet conjugate to the position x'^. Plugging 



(13 1 into (17 1, one straightforwardly finds that 

TTf h — 



dz h\l L^h- 7r| 



(18) 



where tt^ is a constant of integration. In order for ^(z) to be real, the right hand side 



of ( 18 1 must be real. There are three ways this can happen in a manner consistent 
wit 



ith the assumption of steady-state behaviorj^ 



'^Technically, there is a fourth way, but its significance is obscure to us. A string can lead down to 
2 = 2* in the shape of the trailing string, | |20| l, and then turn around and retrace its path back 
up again. The energy localized at the kink must grow linearly with time, which means that this 
configuration is not quite a steady-state solution. 
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One can choose ttj = 0. This leads to ^ = 0, which shows that ([9| is formally a 
solution of the equations of motion. But it is not a physical solution — at least, 
not in the context of classical motions of a string — because of the problem with 
the signature of the worldsheet metric. The action for this solution is complex. 
One can arrange for the worldsheet never to go below z = z*, and choose tt^ 



small enough that the denominator inside the square root in (18 1 is always 
positive. There is indeed a one-parameter family of such solutions, and they 
describe a heavy quark and anti-quark in a color singlet state propagating 
without drag (at the level of the current treatment), one behind the other. 
Similar states were studied by other groups jMlIIlII but they do not capture 
the dynamics of a single quark propagating through the plasma, so we do not 
consider them further here. 
One can choose 

7r^ = ±^yM^-±^=iL=^, (19) 
^* vl — 



so that the denominator inside the square root in ( 18 1 changes sign at the same 
value of z as the numerator, namely z — z^, rendering the ratio inside the 
square root everywhere positive and finite. Choosing the sign that makes ttj 
positive means that the string trails out in front of the quark instead of behind 
it. Although this is technically a solution, it does not describe energy loss and 
should be discarded. Choosing tt^ negative leads to the trailing string solution 



that we are interested in. Eq. (18 1 then straightforwardly leads to 



where we have introduced a rescaled depth variable. 



z 



y=—- (21) 

zh 

There is another way to justify the choice of the minus sign, not the plus sign, in 
(ITglPThe solution (pOl IS non-singular at the future event horizon, whereas the 



solution one would get with the opposite sign choice is singular. To understand 
this point, it is convenient to pass to Kruskal coordinates, defined implicitly by the 
equations 

UV = -llJ^g-^*™"' y _e''*/^« . (22) 

In the region outside the horizon, U <Q and y > 0, while in the region inside the 



future horizon, U > Q and > 0. The trailing string solution (13 1, with given 



by ( 20 ) , can be extended to a non-singular solution over the union of these two 



regions: 



Xi = ^\ogV + vta.Tr^y . (23) 
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Thus, the logarithmic singularity in ^ at ?/ = 1 (meaning z = zh) is a singularity 
not at the future horizon, which is at C/ = 0, but at the past horizon, which is at 
V — 0. Reversing the sign choice in (19 1 would lead to a solution that is singular 
at the future horizon but not the past horizon. Causal dynamics in the presence 
of a black hole horizon can generally be described in terms of functions which are 
smooth at the future horizon. 

At the level of our presentation, it has been assumed rather than demonstrated 
that the trailing string is a stable, steady-state configuration representing the late 
time behavior of a string attached to a moving quark on the boundary of AdSs- 
Schwarzschild. In fact, this has been fairly well checked 

The description we have given of the trailing string is not limited to 
Schwarzschild black holes in AdSs. One may extend this analysis to various other 
black hole geometries which asymptote to AdSs near their boundary) ^^ l ^^ l ^^ l ^^ l ^^ l ^° l 
These geometries describe theories which are deformations of SYM. The literature 
also includes a discussion of the distribution of energy along the string^ ^ and an 
interpretation of the shape of the string in terms of a rapid cascade of strongly 
coupled partons.l^ 



4. The magnitude of the drag force 

As we explained in the previous section, an infinitely massive, fundamentally 
charged quark moving at speed v in the xi direction through an infinite, static, 
thermal medium of A/" = 4 super- Yang-Mills theory can be described at strong 



coupling in terms of the trailing string solution ( 20 1 . The quark cannot slow down 
because it has infinite mass. However, it does lose energy and momentum at a finite, 
calculable rate. In five-dimensional terms, this energy can be thought of as fiowing 
down the string toward the black hole horizon. In four-dimensional terms, energy 
and momentum emanates from the quark and eventually thermalizes. To calculate 
the four-momentum Apm delivered from the quark to the bath over a time At, one 
can integrate the conserved worldsheet current of spacetime energy- momentum 
over an appropriate line-segment 2 on the worldsheet. 2 should cover a time interval 
At, and it can be chosen to lie at a definite depth zq in AdSs. The four- momentum 
Apm i^0 

Apm = - J dty/^p""^ . (24) 

Because the trailing string is a steady-state configuration, four-momentum is lost 
at a constant rate: 

dpm 



-^9P'n.. (25) 

"^There is an explicit minus sign in |24jl which doesn't appear in the analogous equation of one of 
the original works.^l This is due to use of the z variable, which increases as one goes deeper into 
AdSs, instead of the r = L/z variable, which increases as one goes out toward the boundary. 
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In particular, the drag force can be defined as 

dt ^ ' 2TTzla' 



^d.a, = ^ = -^gp\ = . (26) 



Using te| and (151, one obtains 



^drag-^^T^^^. (27) 

If, instead of an infinitely massive quark, we consider a quark with finite but large 
mass m, then using the standard relativistic expression 



leads to 



Fh -_^t2^ (291 



It has been explainecP^ that ( 28 1 receives corrections when a finite mass quark is 



described as a string ending at a definite depth z = 2:* on a D7-brane. While these 
corrections are interesting, it would take us too far from the main purpose of this 
review to give a proper explanation of how the D7-branes modify the physics. Our 



discussion is formal because we derived the result (27) in the strict m = 00 limit 
and then applied it to finite mass quarks. 



From ( 29 1 it is clear that the drag force causes the momentum of a quark to fall 



off exponentially: 

p(t) = ^(0)6-*/*'""- where t^^,,^^ (30) 

In order to make a physical prediction for QCD, we must plug in sensible values for 
m. A, and T. The effective quark mass in the thermal medium is already non-trivial 
to specify precisely, but rric = 1.5 GeV for charm and mf, = 4.8 GeV for bottom 
are reasonably representative values which were used in a recent phenomenologi- 
cal study.l^ We will review here two approache^^U to specifying A and T. The first 
approach, used earlier in a calculation of the jet-quenching parameter q from a light- 
like Wilson loop in A/" = 4 super- Yang- Mills theory,'^^ is to identify the temperature 
TsYM with the temperature Tqcd, and then identifying the gauge coupling ^ym 
of super- Yang-Mills with the gauge coupling of QCD evaluated at temperatures 
typical of RHIC. Because of the proximity of the confinement transition, has 
substantial uncertainty. A standard choice is as = 0.5, corresponding to gs ~ 2.5. 
With the number of colors A'' set equal to 3, one finds A « 67r. We will refer to this 
as the "obvious scheme." 

The second approach, called the "alternative scheme," is based on two ideas. The 
first idea is that it may make more sense to compare SYM to QCD at fixed energy 
density than fixed temperature. SYM has e cx T**, and so does QCD, approximately: 
this approximation is surprisingly good for T > 1.2rc, according to lattice data.l^ 
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Table 1. The obvious and alternative schemes for comparing SYM and QCD. 



Scheme 


Tqcd 


TsYM 


eqCD 


eSYM 


A 


rUc 




tc 


h 




MeV 


MeV 


GeV/fm^ 




GeV 


GeV 


fm 


fm 




obvious 


250" 


250 


5.6'' 


15 


Qiv 


1.5 


4.8 


0.69 


2.2 


alternative 


250" 


195 


5.6'' 


5.6 


5.5 


1.5 


4.8 


2.1 


6.8 



''We set Tqcd = 250 MeV because this is a typical temperature scale for heavy ion collisions 
at v^slvi^ = 200 GeV. 

''We use e/T* 11 for QCD. 



But the constant of proportionality is about 2.7 times bigger for SYM than for 
QCDQ That is, SYM has about 2.7 times as many degrees of freedom as QCD 
above the confinement transition. So esym — eqcd implies Tsym ~ ^000/(2. 7)^/''^. 
This identification leads to a suppression of Fdrag by a factor of v^2?7 relative to 
the obvious scheme. The second idea behind the alternative scheme is that the 
't Hooft parameter A in string theory can be determined by comparing the static 
force between a quark and an anti-quark, as calculated in string theory, to the same 
force calculated in lattice gauge theory. The string theory calculation is based on a 
U-shaped string connecting the quark and the anti-quark. This string pulls on the 
static quarks in a fashion that is similar to how the trailing string pulls on a moving 
quark. The lattice calculation is based on computing the excess free energy due to 
the presence of an external quark and anti-quark in a thermal bath. There is a 
significant difficulty: in the simplest string theory calculation, based only on the U- 
shaped string, the force between the quark and anti-quark vanishes for separations 
larger than some limiting distance r*, and this distance is quite small: r.^, « 0.24 fm 
when Tsym ~ 195 MeV (corresponding to Tqcd = 250 MeV). It has been pointed 
ouiP^ that exchange of closed strings between two long strings describing the quark 
and anti-quark at separations r > contribute to the quark-anti-quark force at 
the same order in iV as the U-shaped string. Unfortunately, it is hard to compute 
the contribution of closed string exchange. The approach'^- is therefore to match 
the U-shaped string computation to lattice data'^® near the limiting distance r^, . 
The result of this matching is A « 5.5. 

As we show in table [T] heavy quark relaxation times tc and ti, are remarkably 
short when one uses the obvious scheme, and somewhat larger in the alternative 
scheme. The uncertainties in tc and tt are substantial: even if one accepts the ideas 
behind the alternative scheme, one should probably regard the resulting relaxation 
times as uncertain by a factor of 1.5. 

An experimental study ■^^ favors a model^ in which tc is roughly 4.5 fm at 

"This mismatch has previously been statecfi3]as a factor of 3 rather than 2.7. Some uncertainty ex- 
ists on both the SYM and the QCD sides, because of finite coupling effects and time discretization, 
respectively; but 2.7 is probably closer to the true figure. 
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Tqcb = 250 MeV, as estimated from plots from a detailed exposition of that 
modelPThis seems to indicate that the string theory estimates of and tb, even in 
the alternative scheme, are too short. However, the results of a recent phenomeno- 
logical stud}^ favor a range of parameters that is consistent with the string theory 
predictions translated to QCD using the alternative scheme. 

Let us briefly review the recent study.l^Sl The starting point is the Langevin 
equation, which in the Ito discretization scheme takes the form 

Ax{t) = ^At Ap{t) = -TpAt + ^(t) . (31) 

Here ^(i) is a stochastic force, assumed to be Gaussian and uncorrelated from 
one time-step to the next. The strength of the stochastic force is related to the 
drag coefficient T by demanding that the relativistic Maxwell-Boltzmann distri- 



bution is preserved by the time evolution (311. Ordinary relativistic kinematics, 
E = yj-p^ -|- TO^, are assumed. It is also assumed that 

r - 7— , (32) 
m 

where 7 is a dimensionless quantity with no p dependence. Evidently, tq^ark — 



1/r. The temperature in (32 1 is Tqcd, whereas the temperature in (30 1 is Tsym 



Comparing these two equations, one finds that the string theory prediction is 

_ ■K\f\ /TsymV _ 1 6-8 obvious scheme ^^^^ 
2 V^QCD/ [2.2 alternative scheme. 



The alternative scheme value in ( 33 1 is fractionally larger than the one quoted in 
the study under discussion ,^21 due to the use here of the factor 2.7 for the ratio of 
degrees of freedom between SYM and QCD, as compared to 3 in previous work.l^^l 
(It is a numerical coincidence that the dimensionless factor 7 is the same, in the 
alternative scheme, as tc in femtometers when Tqcd = 250 MeV.) 

The next step of the studj^^H is to compare Langevin dynamics of heavy quarks 
in a hydrodynamically expanding plasma to PHENIXpSl and STAIft^^ data on the 
nuclear modification factor Raa for non-photonic electrons — meaning electrons and 
positrons coming from decays of heavy-quark mesons. Because of the treatment 
of hadronization, the theoretical results are deemed trustworthy only when the 
transverse momentum pT of the non-photonic electron is at least 3 GeV|^For fairly 
central collisions (impact parameter h = 3.1 fm), agreement between theory (using 
the alternative scheme) and experiment is best for 7 between 1 and 3: see figure [l] 
Thus the prediction ( |33[ ) of string theory in the alternative scheme can reasonably 
be said to agree with data to within the uncertainties of the calculations. These 
uncertainties stem in large part from the difficulty of comparing SYM to QCD; 
however, it is also clear that the treatment of hadronization is a significant hurdle.'^^l 

'The electron carries only a fraction of the of the charm quark that led to its production. This 
fraction varies, but a reasonable rule-of-thumb value is 1/2. 
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(a) b=3.1fm, c+b->e* Z]l°i 

-Y=3.0 

oPHENIX(0-10%) 
• STAR(0-5%) 




"0 123456789 10 
[GeV] 



Fig. 1. (Color online.) Comparisoiffl between experimental data for gold-gold collisions with im- 
pact parameter b = 3.1 fm and theoretical predictions based on Langevin dynamics of heavy quarks 
in a hydrodynamically expanding plasma. The open circles correspond to central collisions in a 
PHENIX experiment, El] while the solid circles come from a STAR experiment.@2! Theoretical pre- 
dictions are plotted for 7 = 0.3 (red), 1 (blue), and 3 (green), with 7 is as defined in l |32| . Reprinted 
with permission from the authors from Y. Akamatsu, T. Hatsuda, and T. Hirano, "Heavy Quark 
Diffusion with Relativistic Langevin Dynamics in the Quark-Gluon Fluid," 0809.1499. 



Fluctuations of the trailing strinJMMHEll provide direct access to the stochas- 
tic forces in (31). In the non-relativistic limit, the size of these forces, relative to 
the drag force, is exactly what is needed to equilibrate to a thermal distribution. 
Indeed, the original calculations of drag forc£3Illand stochastic forced were done 
independently. For relativistic quarks, the stochastic forces are enhanced by pow- 
ers of 1/(1 — v'^), including enhancement of longitudinal stochastic forces by the 
startlingly large factor 1/(1 — u^)^/^. This is larger by 1/(1 — w^)'^/'* than what is 
needed to equilibrate to a thermal distribution. Another issue is that the correlation 
time tcor for these stochastic forces grows with velocity: based on results for the 
relevant Green's function^'' one may estimate 

~ ^T(i_„2)i/4 • (34) 

If we use tcor < iquark ^s a criterion of validity for the Langevin dynamics, then we 
get a limit on the Lorentz boost factor: 

<X^- (35) 



Approximately the same inequality arises from demanding that when the string 
ends on a D7-brane, its endpoint should not move superluminally. The same in- 
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equality, but with m replaced by some fixed scale /i, also arises from demanding 
that the worldsheet horizodHUHl should be at a depth in AdSs corresponding to 
a scale fi where the dynamics of QCD is far from weakly coupled. Plugging in 
numbers in the alternative scheme, with fj, = 1.2 GeV and Tqcd — 250 GeV, one 
finds that l/^/l — v"^ < 30. This inequality should be understood as quite a rough 
estimate, because of the quadratic dependence on the quantity fi which is only 
qualitatively defined. For a charm quark, the corresponding limit on the transverse 
momentum of the non-photonic electron is pr ^ 20 GeV. This is well in excess of 
the highest momentum for which there are statistically significant data; moreover, 
for Pt more than a few GeV, bottom quarks dominate, and for them the bound 
l/Vl — ^ 30 translates to pT 70 GeV for non-photonic electrons. The upshot 
is that heavy quarks at RHIC do not obviously fall outside the regime of valid- 
ity of a self-consistent Langevin treatment based on the trailing string — except for 
the troublesome scaling of longitudinal stochastic forces as 1/(1 — u^)^/*, whose 
consequences, we feel, are ill-understood. 

In another study,'^it is argued that energy loss in the perturbative and strongly 
coupled regimes have experimentally distinguishable signatures for large enough 
transverse momenta, px ^ which will be attained at the LHC. A convenient 
observable that distinguishes between predictions of string theory and perturbative 
QCD is the ratio of the nuclear modification factors for b and c quarks. 



In particular, the drag force formula (29 1 implies that at large enough p^, one has 



i?l'ds- — -0.3. (37) 



where we used the bottom and charm masses quoted right after equation (30 1. In 
contrast, perturbative QCD predicts that 

<QCD - 1 - (38) 
Pt 

at large pr, where Pct is a relevant momentum scale. So according to perturbative 
QCD, i?^** should approach unity at large px- At sufficiently large px, the trailing 
string treatment presumably fails, and perturbative QCD presumably is correct. 



But as discussed following ( 35 1 , it is difficult to give a good estimate of the char- 
acteristic value of Pt where the trailing string fails. Absent a reliable estimate of 
the characteristic pt, the upshot is that if the measured i?^'' is significantly below 
the perturbative prediction, the trailing string should be considered as a candidate 
explanation. 

5. The perturbed Einstein equations 

Given that an external quark dual to the trailing string described in section |3] 
experiences drag, one might ask what happens to the energy that the quark deposits 
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Fig. 2. (Color online.) A visual summarjES! of the calculation of the response (Tmn) of a strongly 
coupled thermal medium to a heavy quark moving with a speed v via AdS/CFT. hmn is a 
perturbation of the metric caused by the trailing string. 

in the medium. At scales much larger than the inverse temperature, one expects 
the excitations present in the medium due to interactions with the moving quark to 
be well-described by linearized hydrodynamics. Earlier investigations of linearized 
hydrodynamics revealed the presence, for generic sources, of both a sonic boom 
and a diffusion wake.l^ The sonic boom is a directional structure, which, in an 
ideal fluid, is concentrated on the Mach cone, but in a real fluid there is broadening 
because of viscous effects. It appears only when the probe is moving faster than the 
speed of sound in the medium, and it comes from constructive interference among 
spherical waves sourced by the quark along its trajectory. The diffusion wake is a 
flow of the medium behind the quark in the direction of the quark's motion. It too 
is broadened by viscous effects. We will discuss the hydrodynamic limit in more 
detail in section [6.1.21 

Apart from a qualitative understanding of energy loss at large distances, little 
is known a priori about what happens, for example, at small distances close to 
the quark. An all-scales description can be achieved using the gauge-string duality, 
where one computes the disturbances in the stress-energy tensor due to the presence 
of the quark. This was done in a series of papers E ^ l ^^ l '^^ l ^^ l ^^ l ^^ l ^^ l ^'^ The linearized 
response of the lagrangian density in the dual field theory was also computed in 
Fourier space IHEH The purpose of this section and the next two is to present a 
reasonably self-contained summary of how the gauge theory stress-energy tensor is 
computed starting from the trailing string. A visual summary of the main elements 
of the computation is shown in figure [2j 

In the context of the AdS/CFT duality, the stress-energy tensor {Tmn) of the 
boundary theory is dual to fluctuations of the metric in the bulk. So in order 
to compute the expectation value of (Tmn), one first needs to compute to linear 
order the backreaction of the string describing the quark on the metric. Non-linear 
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corrections to the Einstein equations will be suppressed by VX/N. The total action 
describing both the string and the metric is S" = 5'buik + *S'string, which one can write 
as 

/^(i? + 12/L2) _ 1 
2^2 2na' 



S = 



d^3 



(39) 



where x'i{a) is the embedding function of the string in AdS5-Schwarzschild. In a 
gauge where we parameterize the string worldsheet by a" — (t,z), x'i{a) is given 

by 



xjf ^ {t vt + S,{z) z) 



(40) 



with ^(z) as given in (20 1. The equations of motion following from (39) are just 
Einstein's equations: 



2 



(41) 



where 



2™' L3 



-^yWl-v^Six' -vt- e(z)) S{x^)6ix')da.x':d''x: (42) 



is the bulk stress-energy tensor of the trailing string. 

To compute the backreaction of the string on the metric, we write 

G — (7^") + h 



(43) 



where g|S is the unperturbed AdSs-Schwarzschild metric given in (|3|, and plug this 
into (41 1 to obtain the linearized equations of motion for the metric perturbations 
h^i,. The resulting equations take the form 



where h^'^ — G'^Q^GI^s^h per and C is the differential operator given bjJ 

Ch^''' = -Dh^, - 2R^p,,hP'' + 2RP^K)p - V^Vh + 2V^P-Vph''^P 

^ ^ R\ hP" 



+ G^o") (-VpV./i"" + Uh + Rp^hP^) 



6 

Z2 



(44) 



(45) 



The covariant derivatives, the Riemann and Ricci tensors, and the Ricci scalar 



appearing in (|45| are computed using the background metric G^^}, and we have 
denoted h = G'0h^"'. 



For a steady-state solution of (44 1, hp^ depends on x^ and t only through the 
- vt. We therefore pass to co-moving Fourier space variables by 



combination 
writing 



|^Tj^''(z)e4^i(^'-''*)+^^^'+^=^"1/^« , (46) 
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where we have defined the dimensionless wavevector K = kzn = k/nT. We make 



a similar expansion for /i^^. In Fourier space, (44) can then be written as 



K 



(47) 



with t'^ being given by 



'K 



lira' Vl - 



V 





V 












h 


h 


V 

h 










































h 










v"^ — hj 



(48) 



where y is the rescaled depth coordinate defined in (21 1, but the tensor components 
are given in the {t, , x'^ , x^, z) coordinate system. The exphcit form of is too 
compHcated to be reproduced here. 

We will decouple equations (47 1 by passing to a gauge where /i^^ 0, which 
we will refer to as "axial gauge." Note that h^z = leaves some residual gauge 
freedom. We will discuss this shortly. 

The rest of this section is organized as follows. In section [STT] we explain how to 



decouple equations (47). In section 5.2 we explain the boundary conditions needed 
to solve these equations. Lastly, in section [5T3] we explain how the one-point function 
of the SYM stress-energy tensor is related to the asymptotic behavior of the metric 
perturbations near y = 0. 



5.1. Metric perturbations in axial gauge 



= — t 



e where h^z 


= 0, 


and let 


/Hqq 




H02 


Ho3 


o\ 


Hio 


Hii 


H12 


Hi3 





H20 


H21 


H22 


H23 





H3Q 


H31 


H32 


H33 





V 











0/ 



(49) 



Rotational symmetry around the direction of motion of the quark allows us to set 
K = {Ki K±_ 0) with if_L > 0. Defining 



K=JKi 



Ki 



d = tan 



(50) 
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one can form the following linear combinations of metric perturbations: 



A- 

Bi 
C 
Di 

El : 
S4: 



-Hii + 2 cot dHi2 - cot^ dH22 + csc2 §H: 



33 



V 



B2 = - 



2^2 

i/i3 + tan i?i723 



— sin i9iJi3 + cos ^H23 
Hqi - cot ■dHQ2 
2^ 



Do = 



-ffii+2cot 2^Hi2 + H22 



2^;2 



1/3 

2 ( ^-^Hqo + Hii + _ff22 + ^33 
^^11 + H22 + H33 



Hqi + tan -dH, 



02 



2w 



-Hii -H22 + 3 cos 2z9(-gii + H22) + 2^33 - 6 sin 2^Hi2 

4 



(51a) 

(51b) 
(51c) 
(51d) 

(51e) 

(51f) 

(51g) 



Using these new variables, the Einstein equations (47 1 decouple into five sets 



.mi 



A = ^e-*-^i«/^« (52a) 
h 



\ K^f-h v^hcosH 



B2 



(52b) 



(52c) 



3 /i' 



9l+[~7 + T)9y + j;2i^'^os'^-h) 



C = 



(52d) 



_3 , V 



f-h v^hcos^d 

/l2 V— 1 ZJ^COS^I? 



/i li 



(52e) 



-hD'2 = Q 



(52f) 
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3 I 3h' 
y ^ 2/1 






V 







/-2h 














V 2h 



-1 + 



12^2 cos^ -d 



-121)2 COS^ I? 



2h 



-3 



6i)2 cos^ 7? 
2/i 
-2/i 




2h 




/i 

-ft 
cos^ i) + hj 









y 




E2 






E3 









( 



y_ -iKti/zH 
h 



1+3 cos 2^9 2 
1 c/ 



(52g) 











-h -3w2 cos^ z9 - /i -h dy 



h 







1 

6h 



-6/i' -3ft' 

-3ftft' ISw^cos^i^ft' S{3v^ cos^ t} + h)h' 
^2K^yh ~2K^v^vcos^d -2K'^y{3v^ cos^ i9 - h) 2K^yhj 









E2 




E3 




\eJ 



ft' 



—ivy seed 

.r. -e-*-^i«/^« I 3ii;ycosi?(w2 + ft) | . (52h) 
^^y^ ^ K{v- - ft) 



In (52 1, A, Bi, C, Di, and Ei are all functions of the rescaled AdS depth y = z/zh, 
and primes denote derivatives with respect to y. 

A few comments are in order. Let's start by counting the equations. Einstein's 



equations ( 47 1 consist of fifteen linearly independent equations that split between 
the A, B, C, D, and E sets as follows: the A and C sets each consist of one second 
order equation for one unknown function; the B and D sets each consist of two 
second order equations and one first order constraint for two unknown functions; 
lastly, the E set consists of four second order equations and three first order con- 
straints for four unknown functions. At first glance, the B, D, and E systems of 
equations might seem over determined. A more careful analysis shows that the con- 
straints are consistent with the second order equations in the sense that if they hold 
at a particular value of y, they continue to hold at all y. We can therefore think of 
the constraint equations as reducing by five the number of integration constants in 
the second order equations. There are fifteen remaining integration constants that 
are fixed by the boundary conditions which we discuss in section [5. 2[ 



Because the Bi and C equations, (52b)-(52dl, have no source terms, one can 
consistently set Bi = B2 



C = 0. This identification is in fact enforced if we insist 
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that the response of the medium should respect the same axial symmetry around 
the direction of motion of the quark that the trailing string does. In order to keep 
the discussion of perturbations, boundary conditions, and integration constants 
general, let us not discard the Bi and C fields just yet. A more general source 
would force them to be non-zero. 

Equati ons (|52[ ) can be reduced to just five equations using the residual gauge 
symmetryl^Il^ The action (39 1 is reparameterization invariant, and infinitesimal 
gauge transformations act by sending 



hn 



(53) 



for any one-form C,^dx^^ . This allowed us to pass to axial gauge in the first place: 
if we had started with some arbitrary metric perturbations h^v we could solve, at 
least locally. 



(54) 



so applying the infinitesimal gauge transformation (53 1 gives h^^ — 0- Residual 
gauge symmetry arises from the fact that ( 54 1 specifies only up to five integration 
constants, which in general are functions of t and x. Put differently, there are five 
linearly independent gauge transformations that preserve the axial gauge condition 
hp,z = 0, and the corresponding are given by the linearly independent solutions 
to (|54| with = 0. 



For the steady-state solution of ( 47 1 the allowed gauge transformations are the 
ones where the t and x dependence of is of the form e*^'"^'"/^*^ , where we have 
defined 



K,n={-K,v K2 K3 



(55) 



As in previous sections, our convention is that lower case Roman indices m,n,... 
take values in (i The first four of these gauge transformations are pa- 

rameterized by 



Ha) 



zne 



27ra' 



where a = 0, 1, 2, or 3, while the fifth is given by 



c, 



^2 ^^^tK^x^/zH 



(5) 



2yVh6'^ - iK 



oV 



2T:a' LV I — V 
The corresponding pure gauge solutions are given by 



(5q — iK^ arcsiny 



(56) 



(57) 



a = 0,1,2,3 

Hmn(5) = 4:Vh7],nn - 2KmKn arcsiu - iy^VhS^^S^ 



2K°K, 



(mVn)0 



arcsm y 



y 



(58) 



(59) 



where r/mn is the Minkowski metric. One can straightforwardly check that (58h(59) 
satisfy equations (47). 
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Using the definitions (51 1, we can work out the pure gauge solutions in the 
ABODE variables. We find that A and C are invariant, which was to be expected 
since their equations of motion are already fully decoupled. The other variables, 
however, do transform non-trivially under (53 1 with (56 1 and (57 1: for example the 
Bi and Di variables vary by 



5Bi 
SB2 = 



5Di ^ 
5D2 = 



iKi{\i - X2 cotd) 
2v 

iKi{\i - A2 cotz?) 



(60) 



w^cos^i? cos^ ^ sin^ i3 

where Xa are arbitrary constants multiplying the pure gauge solutions parameter- 
ized by C^^y 

From these transformation laws it is easy to see that there is a linear combination 
of Bi and B2 that is gauge-invariant (and the same is true for Di and 02). A choice 
of gauge invariants is given by 



B = Bi- cos^dv^B2 
D = Di-cos^dv^D2 



(61) 



The transformation law for the Ei variables is simple to derive but its exact form 
is not very enlightening. The corresponding gauge-invariant can be taken to be 

E = 4:E3 + Ei + {I2E2 - 3Ei)v^ cos^ ^ - {2Ei + Ei)h{z) . (62) 



From (52) we obtain the equations of motion for the B, D, and E invariants: 
'3h^ + {h-A)v^cos^'&\ ^ K^iv^cos^-d-h) 



dt 



yh{v'^ cos^ -d — h) 

r + {h ^ A)v^ cos^ d 
\ yh{v^ cos^ d — h) 



B^O 



f 1 - COS^ I? + 



a, 



/l2 

Aivy^ 



D 



K{v^ cos2 'd-h) 



(63) 



(64) 



4 16^^^ 
yh 4 - 6u2 sin^ {) + 2h 
Klv^ 32y« 

" y2/j(4_6u2sin2l? 

-'^le y 31)2 (,Qg2 ^ _ 2 - i;2 



2h) 



h cos2 d 



E 



(65) 



h 2 + h-3v^ cos2 d 
9v^ cos* ^ - 18v^ cos^ 1? 



1 



3vKi 



y^ + 9 



The gauge invariants we just described are unique up to an overall z-dependent 
factor, provided we only consider linear combinations of the Bi, Di, and Ei vari- 
ables. If we also allow derivatives of the latter there are many other choices of gauge 
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invariants that might prove useful. For instance, 



= -K'^E 
" 3 



2v^sm'dzHD[ 
vzhB[ 

-Asin^ 

-C 



At 



(66) 



is another set of gauge invariants which has been used in the hteratureP^l Yet 
another set, consisting of "master fields,"'^ is: 



2 

-^J— — hD^ sm tf 



[3h^yE2 + 6h^E'2 - 3h{3 + y^)E'^ 



2hyK^ {El - E3) + 3y (2%^ _ k^) {AE2 - E^)] , 



(67) 



where ay — l/27ra'\/l — V'^. The gauge invariants (67) will prove useful for the 



asymptotic analysis of section |6.2[ Below, we list their equations of motion: 



3 h'\ K^v^ cos'' ^ - h) , 
--y+h)'^-+^^^ 



3 h'\ , K'v'y' cos' d + h{3 - K'y^ + 9y-^) 
-y + j)^'^ + y^^ 



Kfv' K'{A- K'y') + 12y'{6-y^) 



h' 



y^h{K^ + 6y2)2 



-Jt (68a) 
-Jv (68b) 



- -Js (68c) 
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where 



Todd 



Tcvcn 



TCVcn 



Todd 



-iKi^/zH 



h 



\y sin i?(5j/'* — 1 — iKvy^ cos "d) + iKy'^ tan t?] 



6i4 (i^2 6y2)2 
3X^(2 - ?)v'^)y^ + 18(2 + 'iv'^)y^ - 3iKv^y^{K^ - 12?/^) cosi? 



+ 9iKv^y^{K^ - ISy^) cos^ + 9ii:^w''y^(i^2 + Gy^) cos" 
-6ivy^K{K'^ -2y^)sec'& . 



(69) 
(70) 

(71) 



(72) 



By omitting the "even" and "odd" superscripts in (68a I and (68b) we mean that 
these equations take the same form (up to the different source terms as indicated 
in ©-((Til). 



5.2. Boundary conditions 



In axial gauge, the system of equations (47 1 consists of ten second order differential 



equations in y and five first order constraints, so we need to specify fifteen integra- 
tion constants. The purpose of this section is to show that five of these are fixed 
by imposing boundary conditions at the horizon (y = 1) and ten of them are fixed 
from the boundary conditions at the conformal boundary (y = 0). 



Consider first the tensor field A, defined in (51a I. The near horizon solution to 



(52a I takes the form 



^(;T-ln4) 



4 1 



2 , 



(1 - y)i-™^^/Vc/A(l-2/)-™'^^/''+K4(l-y)™'^^/' + . 



(73) 



where . . . indicates terms which are subleading to one of the ones shown. The 



first term in ( 73 1 is a particular solution to ( 52a I characterizing the response of 



A to the trailing string source. The second and third terms are solutions to the 
homogeneous equation, and Ua and Va are constants of integration. The correct 
boundary condition at the horizon is Va = 0. This corresponds to requiring that 
there are no outgoing modes at the horizon, and it can be justified by the fact that 
classical horizons don't radiate. The result of choosing purely infalling conditions 
is that in the dual gauge theory, we describe a causal response of the medium to 
the probe. 



To see that the Ua term in (73 1 is infalling, let's define a new coordinate y^ — 
log(l — y) that ranges from — oo when y — I (the horizon) to zero when y — 
(the conformal boundary of AdSs). Recalling that the time dependence of metric 



24 S. Gubser, S. Pufu, F. Rocha, and A. Yarom 



perturbations was assumed to be e ^^^i^^/^h^ Qj^g gggg immediately that the (1 — 
y'^-ivKi/A i^gj-jj^ (73) corresponds to certain metric components behaving as 

/imn~e-™^i(^*+4*/"«)/4 (74) 

at large negative . This behavior describes a wave traveling towards negative y* , 



i.e., falling into the black hole horizon. Similarly, the Va term in (73 1 corresponds 



to an outgoing mode. The same story goes through for the C combination of metric 
components defined in (51c I, except that since there is no source term, there will 



be no analog to the first term in ( 74 1 



A subtlety arises in the horizon boundary conditions for the B, D, and E sets: 
for each set, in addition to a single infalling solution and a single outgoing solution. 



there are the pure gauge solutions discussed around ( 58 1 and ( 59 1 . The pure gauge 
solutions are neither infalling nor outgoing at the horizon. The correct boundary 
conditions are to exclude the outgoing solution and to permit both the infalling 
solution and the pure gauge solutions. If one passes to a description only in terms 
of gauge-invariants, then this subtlety is avoided: each gauge invariant field has 
only an infalling and outgoing solution, and the latter is excluded by the horizon 
boundary conditions. 

Having fixed five integration constants (one for each of the ABODE sets) using 
infalling boundary conditions at the horizon, we now discuss the boundary condi- 
tions at the boundary of AdS. Close to y = 0, Einstein's equations can be solved 
in a series expansion in y. The two homogeneous solutions are 

H^l{y) = R^,,{l + 0{y')) H^Uy) ^ Q,^^ {y^ + 0{y')) (75) 

where Rmn and Qmn are arbitrary constants. The full solution then takes the form 



(2), 



The components of Pmn are given by 

1(1 + 2v^) 



P„ 



(l{2 + v^) 

-2v 



nin 3 

y 



0{y'). 



V 



1(1 -.2) 





id 



The boundary conditions we impose are that 



(76) 



(77) 



(78) 



thus fixing the remaining ten integration constants. Allowing non-zero Rmn would 
correspond to deformations of the gauge theory lagrangian. With this condition 
imposed, the Qmn are related to the expectation value of the stress tensor in the 
gauge theory, as we will see in section |5.3| 



Equation (78 1 can be translated into boundary conditions for the ABODE 



variables. Once the horizon boundary conditions are also taken into account, it 



follows that each set of equations in ( 52 1 is supplemented by just enough boundary 



Energy loss and the gauge-string duality 25 



conditions to uniquely fix a solution. For example, ( 78 ) implies boundary conditions 



on A which fix one constant of integration. Another integration constant is fixed 



by the horizon boundary conditions. Since the underlying equation (51a I is second 



order and linear, the solution is unique. A more complicated example is the E set, 



where the boundary conditions from ( 78 1 fix four integration constants and the 



horizon boundary conditions fix one more. Since the four second order equations 



of motion ( 52g I are subject to three constraints ( 52h I , the number of integration 



constants available is five. Thus one again finds a unique solution. If we pass to 



a description in terms of gauge invariants, then (78 1 fixes a single constant of 
integration for each gauge-invariant field. Because the horizon fixes another constant 
and the underlying equation for the gauge- invariant is always second order, we again 
recover a unique solution. 

5.3. The boundary stress- energy tensor 

The AdS/CFT duality offers a prescriptiodMI for computing the stress tensor of 
the boundary theory from the bulk action USESESMl 

We now explain what S'totai, 5mn and e are. g^n is the metric on the conformal 
boundary of AdSs. After taking the variational derivative, gmn is set equal to the 
metric of the boundary theory, i.e., the Minkowski metric rjmn in our case. S'totai is 
given by 

S'totai — S'bulk + SgH + S'counter ■ (80) 

Here S'bulk is the bulk action ([T]); Sqh is the Gibbons-Hawking boundary term 

5gh = ^ J d^xVc^K^ , (81) 
where S is a co-dimension one surface close to the AdS boundary with outward 



normal n^, induced metric 



G^, = G^, - n^n, , (82) 



and extrinsic curvature tensor 

K^. ^ -G^p Vn. ; (83) 

and Scountcr is an additional boundary term which renders the on-shell action finite 
for geometries which do not induce a trace anomaly in the boundary theory. An 
explicit expression for this term is 

Scountcr ^^J d'x^^ - ^ R^^ , (84) 

with i?^ the Ricci scalar constructed from G^j„. Usually, there are additional terms 



in (84 1 coming from the matter action. Since we are working in the probe limit, we 
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do not need to worry about these extra termsF^The coordinate e in (79 1 specifies a 
hypersurface S(e) which coincides with the boundary of AdSs as we take the e — *■ 
limit. 

Since we are working with a flat boundary metric, the variation of iScounter will 
not contribute to (Tmn), and we can choose S to be a surface of constant z. The 
outward normal form is 

n^dx^ = ^dz , (85) 

zyh 



and then 



2 dGmn .„„x 



Equation ( 79 1 now reads 



L 



2 



(T„„,) = hm (if^„ - K^Gl„) . (87) 



Using (87 1, the expectation value of the stress-energy tensor in the absence of 
the quark is that of a thermal bath 

(T„„)bath = ^N^T^dmg{3, 1, 1, 1} . (88) 

o 

The presence of the quark generates two additional contributions: writing 

(T ) - (T k .V, + /" KT'^ + (T^ ^1 pnKii^'-"t)+K2X^+K3X^]/zH 

yj-mn/ — \J-mn/bath ^ J ^27^)3 L\ mra/div ^ K-'-mn/i ^ i 

(89) 

and using the definitions ( [49| and ( [76| , as well as the boundary condition Rmn — 
and the AdS/CFT identities ([2| and ([l5|, one obtains 

(-^mra)div ^ (^™" ~ VrnnPl ) (90) 

(T,l) = (Q_ - wQ/) . (91) 



Let's first understand the divergent contribution. Plugging (77 1 into (|90|, it is 
not hard to see that in position space (T,„„)(jiv takes the form of a contact term 

(T™„)div = ^u„,uWl-v^S{x^ - vt)S{x^)6{x^) (92) 

where = y==f (1, iT) is the four- velocity of the quark. A divergence of this form 
was to be expected, and it is associated to having an infinitely massive quark. In the 
dual gravity language, the mass of the quark can be identified with the energy of the 
trailing string. But this energy is both IR and UV divergent. The UV divergence 
corresponds to the "bare mass" of the quark and is exactly given by 
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where e is the UV cutoff. Equation ( 92 1 then tafces the form of the stress-energy 



tensor of a particle with mass M moving with velocity v along the x^-direction. It 
can be shown that any string configuration whose endpoint lies on the boundary of 
AdS, whether it is stationary, moving with constant velocity relative to the plasma, 
or accelerating, will generate a divergent contribution of the form ( 92 ) !^ 



The finite contribution to the stress-energy tensor given in (91) can be further 



simplified by using the 55 Einstein equation. Using the series expansion (76 1 with 
R„ 



0, the 55 Einstein equation imposes a tracelessness relation on the Q„ 



Qoo + Qii + Q22 + Q 



33 



0. 



This equation implies further that 



Qn 



(94) 



(95) 



Note that even from (91 1 one can see that {T:^^) is traceless, so (94) does not imply 



tracelessness of the boundary stress-energy tensor. What (94) does is that it allows 



us to write the one-point function of the stress tensor in the simplified form ( 95 1 
The 5m Einstein equations imply four more relations among the Qmn- 



IV 

2 



{v -10 0). 



(96) 



Using ( 95 1 , ( 96 ) shows that the boundary stress-energy tensor fails to be conserved: 



{v -1 0) 



(97) 



2 ^/^^ 

This non-conservation comes from the fact that the quark is prescribed to move at 
constant velocity. The drag force can be interpreted as minus the force exerted by 
the quark on the medium. With this interpretation in mind, we can check explicitly 



that the drag force (27) can be recovered from (97), as follows. Given the stress- 



energy tensor, the external force acting on a region V can be computed from 
d 



dt 



d^x (T"") 



dV 



Jv 



(98) 



The region V is assumed not to depend on time in the asymptotic rest frame of 



the plasma. The first term in ( 98 1 gives the rate of change of energy- momentum 



in this region, while the second term corresponds to the energy-momentum flux 
through the boundary of V. To obtain the last equality we used the divergence 



theorem. Using ( 89 ) and taking the limit where the volume of V goes to infinity, as 



appropriate for computing the total force on the system, one can see that 

V ttTVA 



F'^^izif lim i^™(Tr 



2 



(w 1 0) 



(99) 



is then indeed minus the drag force, as can be easily checked by comparing (99) 
to (|27|. 



28 S. Gubser, S. Pufu, F. Rocha, and A. Yarom 



The prescription ( 95 1 gives the stress tensor in terms of the coef&cients Qmn 

rnn- We Can 



that appear in the near boundary asymptotics of the components of H 
give similar prescriptions in terms of the near boundary asymptotics of the other 
linear combinations of metric perturbations that we discussed in section |5.1| — all 
that we need to do is to relate the Qmn to a near boundary expansion of these 



variables. For the ABODE variables defined in (51 1, the near boundary behavior 
is 

Px 
3 



X 



Rx+'-^y^ + Qxy^ + Oiy'') where X = A, B„C, D,, E, , (100) 
Here, as before, the Qx and Rx are arbitrary constants while the Px are set by 



( 52 1 . The boundary condition R„ 



is equivalent to Rx — 0. 



It is easy to invert the definition (51 1 and use the near boundary expansions 



( 100 1 and (76) to write the Qmn components in terms of the Qx- However, it is 



useful to notice that we need not specify all of the Qx- The reason for this is that 



the constraint equations (52c 
namely 



3i - Qb2 — Qdi - Qd2 — - 



(52fl, and (52hl imply relations among the Qx, 



4vKi 



Qe, + 2Qes = 



Qei — 2(9 B2 



IV 

2K\ 



(l-3l>2cOs2 ^)QE^+2QE, 



3iv{l + v'^)cos-d 
2K 



(101) 



Using (101 1, we can write Qmn in terms of only five of the Qx, one from each 



set. This will be nothing more than a parameterization of the most general Qr, 



obeying constraints ( 94 1 and ( 96 ) 



Recalling that the response to the trailing string has Bi = B2 = C = every- 
where because of symmetry, we conclude that we can write Qmn hi terms of Qa, 
Qd = Qdi and Qe = Qei |f] 



Qr. 



where 



v"^ sin^ I? 




Id + e„ 







>E+Pn 



0\ 

-2 sin^ d sin 2i9 
sin 2-d -2 cos^ d 



\0 



(102) 



(103) 



2/ 





-2 sin 2?? 



V 



-2 SIT? 2i? 
sin 4z9 




-2sin2i? 0\ 

sin4^? 
2 sin^ 2d 







0/ 



(104) 



sin later sections, we will continue to use Qjj to mean Qo^, ^.nd likewise Qe = Qei- If the 
gauge-invariants D and E defined in ( |61[ l and (j62| were expanded in powers of y, like in ( |100| , 
the coefficients of would be related to Q/j^ and Qei, but not identically equal. 
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-4 


4u cos^ ■& 


2t;sin2t9 





Av cos^ 1? 


4eii 


4ei2 





2v sin 2z9 


4ei2 


4e22 














2^2 cos^ ■& 



eu = ^ (-1 + (1 + v^) cos2 d - cos^ (^O^) 

ei2 = ^ sin 2?9 (l - 3v'^ cos^ ?9) 

622 = ^ cos^ z9 (-1 - 2v^ + 3v^ cos^ d)) 



ycos?? 
Sic 





4w 
4u tan i9 


2\ 



4:V 4w tan ?? 

Pll Pl2 
Pl2 P22 




Pll = (1 - 3w^) cos 2i? - 5 - 
pi2 = (1 - 3w^) sin 2i? - 4 tan d 
P22 = (Bw^ - 1) cos + . 

For the master fields (|67|, the near boundary behavior is 







2 + 2vy 



(106) 



/oven _ 



Pt 3_^Qt 4 
PV 2 



, Qv 3 
2 , 

^2/ 4 



(107a) 
(107b) 
(107c) 



Here we write asymptotics after having imposed the condition Rmn — 0. Similarly 



to before, the Q coefficients are arbitrary while the P are set by (68). The odd 
master fields ip^'^ and V'y^'^ have similar expansions, but since their equations of 
motion are homogeneous, we can set them to be identically zero and worry about 



them no longer. Using the definitions (67) and (101) it is easy to relate the Q 



coefficients for the master fields to the Qx- The result is 



Qd = 
Ob- 



Qt 

a^w^L^z|^ sin^ 



4vKi 
iv 



Qs 



2Ki 2a^L^ZH 



(108a) 
(108b) 
(108c) 



Similar formulas relating the Qx to the asymptotics of the other gauge-invariants 
discussed can be easily derived, but we will not need them. 
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6. Asymptotics 



Equations (52 1 or (68 1 are difficult to solve exactly and we eventually resort to 
numerics to obtain a full solution. However, there are various approximations which 
can be used in order to get a handle on the large and small momentum asymptotics 
of the metric fluctuations. These may be Fourier transformed to real space, giving us 
approximations to the near- and far-field behavior of the boundary theory stress- 
energy tensor. In section 6.1 we focus on the small momentum asymptotics of 
the solution: in section |6.1.1| we construct a small momentum series expansion of 
the metric perturbations, while in subsection |6 . 1 . 2| we explain how this maps into 
hydrodynamic behavior on the boundary theory. In section |6.2| we focus on the 
near- field of the stress-energy tensor. 



6.1. Long distance asymptotics 
6.1.1. Momentum space analysis 

The small momentum asymptotics of the solution can be obtained by formally 
expanding the appropriate fields in power series in the momentum K. In what 
follows we will go over such an expansion in the A, B, C, D, E variables of (51 



A similar construction can be carried out for the other parameterizations of the 
Einstein equation^ZEH given in (66 1 and (67). We start by formally expanding the 



A variable in (51a I such that 



A 



J2 



(109) 



where the a„ are functions of K which are invariant under rescalings K 
That is, the a„ only depend on the direction of K, not its magnitude K 



XK. 

= \K\. 

Plugging the expansion (|109|) into (|52a| and collecting terms with identical powers 



of K, we obtain a set of equations for a„ of the form 



^ y 3 ^n , 
where Sn can depend on a,n with m < n. For example, 

y „ .y cos i3£, y f cos cos^ — /i 
So=T Si = -i-— 02 = -t\ T3 "0- 



The most general solution to (110 1 is 



y y3 ry ^h(§) 



(110) 



(111) 



(112) 



In order to satisfy the boundary conditions ( 78 1 at the asymptotically AdS bound- 
ary, we require that j/o = 0. The other integration constant in (1121 is obtained by 
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matching the near horizon behavior of (1121 to the series approximation (73) with 
Va = 0, expanded at small K. For example, at order n = 0, we find that 



ao 



2 tan + log 



1 + 2/ 



yi log 



i + y' 



1 



(113) 



Expanding ( 113 1 near the horizon and matching it to a small K expansion of (73 1, 

A^il-y) + UA + OiK), (114) 

we find that we need to set j/i = 1. 

According to (95 1 and (1021, to obtain the stress tensor we will eventually 
need Qa- Expanding (1131 for small y and recalling (100 1, we can read off that 
Qa = 1/4 + 0{K). The next order corrections to Qa can be obtained in the same 
manner. We find that to order 0{K), 

1 i In 2 

Qa ~ 1 

^^4 8 

A similar treatment gives us the asymptotic values of Di and Ei. The result is 
i sec "d sec^ •§ — Av"^ 



(115) 



Qd = - 



16u2 



Qi 



3ii;(l + v^) cos 1? 3i;2 cos^ ■d{2 + ^^(l - 3 cos^ i?)) 



(116) 

2if(l-3i;2cos2i?) 2(1 -3^;2cos2|?)2 • ^^^''^ 

Note that Qr, exhibits a pole structure at iiT = 0, while Qe exhibits a pole structure 
at = 3Kfv^. We will see shortly that these correspond to the diffusion pole and 
sound pole expected of the hydrodynamic behavior of the plasma far from the 
moving quark. 



6.1.2. Relating the large momentum asymptotics to hydrodynamics 

At scales much larger than the mean free path, we expect to be able to describe 
a thermal gauge theory by effective, hydrodynamic, slowly varying degrees of free- 
dom. Let's first see what such a description entails, and then check how well our 
asymptotics match it. Consider a conformal theory for which the hydrodynamic 
energy-momentum tensor {T^^'^'^°)mn is traceless: 

^T^ydroyn^Q^ (118) 

A static configuration of the fluid will be given by 

(r''^''™)™„ = |diag{3,l,l,l} . (119) 

If we now perturb this fluid slightly, and choose our hydrodynamic variables to 
be e = (rhydro-joo ^ (yhydro-jOz ^j^j^ i ^ 1, 2, Or 3, then we should be 

able to write the remaining space-space components of the energy- momentum ten- 
sor, {T^^'^'^°)ij, in terms of gradients of e and Si. The only possible combination 



consistent with the tracelessness condition (1181 is 



(yhydro)^^. ^ „ ^Uk^.S,) + 0{k') , (120) 
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where the free parameter T is the sound attenuation length, related to the shear 
viscosity, temperature, and entropy density through T ~ Arj/isT . Using Q we find 
that r = 1 /SttT in the theories we are considering. Brackets denote the symmetric 
traceless combination, i.e.. 



(121) 



The energy density e and the energy flux Si can now be computed from the con- 
servation equations. 



(122) 



where f^'^™ is the source which perturbs the hydrodynamic stress-energy tensor. 
Separating the spatial part of the source, f^y'^''° ^ into longitudinal and transverse 



components, 



/ 



/^^^™, with {k). 



we find from (120 1 and 0221 that 



S-k 
S- {S-k)k 



-3i (/>dro.fc_ jhydro^ 



s/ofc^r 



k^ — 3aj^ — SiTk^Lo 

^f^y'^'^e + s/^'i™ • kuj 



Oifk) 



fc2 — 3uj^ — SiTk^u) 
0{fk). 



■ Fliydro 

'■Jt 



0{fk) 



(123a) 
(123b) 
(123c) 



In real space, the pole at approximately lj Q corresponds to a diffusive wake 
behind the source, and the pole at roughly fc^ ^ iu)'^ corresponds to the shock wave 
which appears along the Mach cone behind the moving probe. The displacement of 
the two poles from the real axis due to the shear viscosity of the fluid is responsible 
for viscous broadening of the wake and Mach cone. 

to the hydro- 



We can now compare the small momentum results of section 6.1.1 



dynamic behavior (1231. To do so, we use (115)-(117l in (1021 and (95) to obtain 



the small momentum stress-energy tensor. We find 



\ mn/ 



(7rT)VA 



e 

\ 



^Si -Si_ 0\ 



7}. 



+ + o{K) 



(124) 
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with 



Si 



27r {K^ - 3Kfv^) 
iKi (t;^ + 1) 

i Kl 



S< =- 



2KiTT mKliTv^ 



2tt {K^-3Kfv^y 
-6Kfv^ - K^v + KlK^ (6w2 + l)v 

K1K2V {K^ + 3K^v^) 
27r(if 2 - SK^v^y 



(125a) 



and 



Tij — gC^jj 



6 



iKiiS 



J. rn 



diag {0,-2, 1,1} 



Note that we may always carry out a resummation 



(where are constants) or 

1 0{K^) _ 1 + 0{K) 



(125b) 

(125c) 

(126) 
(127) 

(128) 
(129) 



After such a resummation, the expressions for e and S in ( 125 1 take the form ( 123 ) 
with 

^hydro^^(_^2 „ 0)^-*X"r„™ (130) 

and r = I/SttT. Thus, if we identify the first term in the parenthesis on the right 
hand side of (1241 with the hydrodynamic contribution to the stress tensor of the 
SYM theory, {T'^^'^''°^^^^)r,in, then we find that (Ti^y^^-SYM)^^ satisfies 

.^m(j,hydro,SYM)^^ = /hydro,SYM (^3^) 



This should be compared to the full conservation law ( 99 1 , 



(132) 



The extra term 7^„ in ( 124 1 holds information on the deviation of the stress-energy 
tensor from its hydrodynamic form. Alternately, J„ — ifc'"7^„ gives us an effective 
hydrodynamic four-force which sources the hydrodynamic stress-energy tensor. It 
is no coincidence that the large distance asymptotics of the stress-energy tensor 
agree with a hydrodynamic expansion. In fact, it can be shown that generic probe- 
sources excite the metric in such a way that the resulting large distance asymptotics 
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of the boundary theory stress tensor will have hydrodynamic behavior .'22' There is 
also mounting evidence that such a connection between hydrodynamics and gravity 
goes beyond the linearized approximationJ^^ 



To see that the pole structures in (1281 and in (1291 really correspond to a 



laminar wake and a shock wave, we Fourier transform them to real space. Consider 



the last two terms in (125b), resummed as in (129|). Using 
d'^K e'^-^ 2 /X 



7i-d/2 



(2^)^ (X2 + /,2)„ {iny/^T{n) \2fi J ^n-dMf^X) , 

we can Fourier transform the resummed expression to position space: 
d^K 2v 



(133) 



(27r)3 7r(if2 _ AiKiv) 



j(Ki{x^-vt)+K2X^+K3X^)/zH 



VZH 



2^2 



vt] 



(134) 



As expected of a diffusion wake, we find that the configuration (1341 exhibits a 



directional energy fiow, with a parabolic shape far behind the moving quark. 



Fourier transforming the resummed sound pole ( 128 1 is difficult due to the cubic 



terms in the denominator, and we eventually resort to numerics to convert such 
expressions to real space. To see that the pole at ^ 3Kfv^ really corresponds to a 
shock wave, it is sufficient to Fourier transform only the leading order contributions 
to these poles. Since we are neglecting the viscous contribution to the pole structure, 
this corresponds to the inviscid limit. Using contour integration and the identitied^ll 



OO Q 

3 o{ax) sm{bx)dx = < 
1 



/■OO 

/ 3o{ax)e~'^^dx 
Jo 



<b < a 
= < a < & 

Re(a ±ib) >0, 



(135) 



(136) 



where J is a Bessel function of the first kind, we find that, for example, the leading 



contribution to the energy density (125a I reads 



d^K 3iKiv{l + V^) ^(^K^{x^-vt)+K2X^+K■,x^)/zH 

(27r)3 27r(i^2 _ 37^:2^2) ^ 

3D(-u^ + l)a;^2|- 

fiTT^[(x^Y + (l-d,V^)x\Y'^ 



f2 < 1/3 



3i;(-u^ + l)a;^z^ 



i-n'^[{x^Y + (l-3v^)x\) 



~SJ2. 







f2 > 1/3, 
otherwise . 



V3w2 _ y < a;i < (137) 



In obtaining (125a I we assumed that the sound poles are slightly below the real 



K\ axis for > 1/3. This assumption is motivated by the fact that the viscous 
corrections that appear at the next order do shift the poles to the lower-half complex 
Ki plane. As expected, since we have treated viscous contributions as infinitesimal 
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in the leading order result (1371, the real space expression for the energy density is 



singular along the Mach cone. 

6.2. Short distance asymptotics 

As was the case for the large distance asymptotics, the short distance asymptotics 
of the solution can be obtained by formally expanding all variables in large mo- 
menta. Starting from the z coordinate system in ([3]), we look at momenta which 
are much larger than the inverse temperature scale, K = kzn ^ 1. In this case, 
it is more practical to use the dimensionless radial coordinate Z = zk instead of 
y = z/zh- Setting _fC 3> 1 means that we're pushing the black hole horizon off 
to Z — !■ oo, implying that we're nearing the zero temperature limit. To see this 
explicitly, consider the expansion 



oo 



V;^™" = K-^ 2^ i„i^-" (138) 



n=0 



of the tensor modes defined in (67 1, similar to (1091. Plugging the expansion (138) 
into (68a I and collecting terms with similar powers of K, we find that the t„'s 
satisfy 

9|t„ - %dztn -(l- t„ = -(Jt)„ (139) 



where, as before, we need to compute {JT)n order by order in a large K expansion. 
For n = 0, 1, and 2 we find 

_(j^)„ = _£l|!^sin2^Z (140) 
-(Jt)i = (141) 
-(JT)2 = -^^sin2^^^. (142) 



To solve ( 68a I we use the method of Green's functions. The homogeneous version 



of ( 139 1 can be easily solved. The solutions are 



fW(Z) = Z%(^/l-^z) t'^^\Z)^Z^K,Ul-^^z\ (143) 



where I2 and K2 are modified Bessel functions of the first and second kind. In ( 143 ) 



we suppressed the index n because the homogeneous parts of ( 139 ) are the same 



for all n. To simplify the notation we define a — \/\ — kfv'^/k'^. The solution which 
vanishes near the boundary, and therefore does not correspond to a deformation of 
the theory, is t'^' . The other solution, t*^^' , is the only solution which does not diverge 
exponentially in the deep interior of AdS. Clearly, t^^^ is the solution which captures 
the boundary asymptotics we have in mind, and it seems physically reasonable to 
disallow a solution which is exponentially divergent at large Z. A more rigorous 
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approach would be to find a uniform approximation to the two hnearly independent 



solutions of the homogeneous version of (68a) and to show that the solution that 



behaves like t^^"^ near the boundary is not purely infalling at the horizon. This can 
be carried out via a WKB approximation.^^ Thus, the solution to 



dlG{Z, Z') ^G{Z, Z') - f 1 - ] G{Z, Z') = 5{Z Z') 



fc2 



(144) 



IS 



and then 



G 



' {z')-H'^\z')t^^\z) z<z' 

{Z')-H^^\Z')t^'^\Z) Z>Z', 



UZ)= / dZ'GiZ,Z%JT)n{Z') 



(145) 



(146) 



For n — and n — 1 these integrals may be carried out exactly.'^ We find 



toiZ) = -- 
ti{Z) = 

t2{Z) = - 



.3 



L6 



■ sin^ d 



Z^La {aZ) - ZH2 (aZ) 



2a , 
3^' 



(147) 



tfjV^ay . 2 a iKivZ^ 
sm w 



L6 



with L2 a modified Struve function. Recalling (107a I we can read off 



-zjrv'a,, sm^ "dL^ 



y3 2, 



1 

16' 



ivK-i 



3{K^ - Kfv^) 



(148) 



from ( 147 1 . The equation of motion for to coincides with the one that would have 



been obtained starting from a string hanging straight down from the boundary 
of AdS space, boosted to a velocity v in the xi direction. The solution to then 
corresponds to the tensor mode metric perturbation in response to this string and, 
as we will see shortly, it captures the near-field physics of the stress-energy tensor in 
response to a massive quark. At scales much smaller than the mean free path, one 
may effectively ignore the interaction of the quark with the plasma. The function 
t2 corresponds to the first thermal corrections to the near-field of the quark. 

The computation of the large K asymptotics for the vector and scalar modes 
follows in a similar manner. Let 



(149) 



similar to (1381. Expanding (68b I at large K, we find that the w„'s satisfy 

3 



dlu 



z'^n - -^dzVn + (3 - Z'^a'^)v„ 



'{Jv)r 



(150) 
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The first few terms in —{Jv)n are given by 



-{Jv)o 
-(Jy)i =0 

-{Jv)2 = - 



' 3X2^1 



The homogeneous solutions to (1501 are 

and using the Green's function method we find 



(151) 
(152) 

(153) 
(154) 



-Z^(Li(aZ) -li(aZ)) 



3XiX2a4 



which, recalling (107b I, implies 



3fc2 - 4fc2z;2 
3ivkki{k'^ — kfv"^) zjj 



^^ (155) 



(156) 



The details of the computation of the scalar modes can be found elsewhereH33HEni 
The final result is 



>s = a« 



IV 



TT (2 + V^) (fc2 - fc2„2) _ „2(i _ ^2)^2 

12v/fc2 - klv^ 
9 fci(5-llw2) 2w2(l -t>2)/c3 



where 



9^1^ V^l k^~V^kl (F _ fc2^2)2 



+ 0(z^^) 



(157) 



(158) 



With Qt, Qvi and at hand we can use (1081, (102), and (95) to obtain the 
leading large momentum asymptotics of the stress-energy tensor. The momentum 
space expressions for the energy density and Poynting vector are 



(2 + i;2)^/^2 _^2„2 i^2^2(^_^2) 



24 

14 + 7^2 



24^i^2 _ ^2^2 

iKfv^{l - i;2) 



iKiv{llv^ -5) 
187r(if2 _ Kfv^) ^ 24 (i^2 _ ^2^2)3/2 ^ 97r(if2 _ Kfv"^ 

Klv^ (I0t;2 - 1) _ Ktv\l-v^) 

24(if2 _ X2t,2)5/2 " 8(^2 - if2^2)7/2 



(159) 
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9w 



3tt{K^ - Kfv^) ^ 16(i^2 _ if2„2)3/2 + 9(i^2 _ Kfv'^)'^TT 
48(if2 _ i^2^2)5/2 ~ 8(i^2 _ K^y2y/2 



(160) 



and 



tt^T^TA KiK2v{l - v^) 
iKfK2V^{l - 1)2) KiK2v{l + Uv^) 



KfK2V^{l 



9-k{K-^ - Klv^Y 48(if2 _ X2„2)5/2 8(^2 - i^2^2)7/2 



(161) 



The expression for (TJ^) can be obtained from (1611 by exchanging K2 with if; 



In real space, using the notation in (91 1, we find that to leading order in xT 



(T 



K 



V J 7 



(162) 



where Am„ represents a Lorentz transformation with boost parameter v in the xi 
direction and T^""'*^ is given by 
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(163) 



^mrf' is the stress-energy tensor of a stationary heavy quark. Up to the overall 
multiplicative factor, it can be determined by the requirement that it is conserved 
and satisfies conformal symmetry. Thus, the leading short distance behavior of the 
near field of our quark is a boosted version of the stress-energy tensor of a stationary 
quark. At distances much shorter than the typical length scale of the fluid, the quark 
does not see the plasma it is moving through, and behaves as if it were in vacuum. 
Of more interest are the subleading corrections to the stress-energy tensor. These 
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are given by 
(Ttt) = 



^/XT^ v{x - vt) [x\{-b + 131)2 - ^v'^) + (-5 + llv'^){x - vtf 



(164a) 



/AT2 v'^{x-vt)[{l-v'^)x\+2{x-vt) 



^/^^ 24[xi(l-i;2) + (a;-i;t)2]5/2 ^^^^^^ 

VAT^ :c^(1 - z;^)^;^ [8xi(l - z;^) + ll(x ~ vtf] 

{-i-tx,)^ . rM (164c) 

Vl-v^ 72[a;2_(l-w2) + (x-vt)2]^/2 

^_ ^/A^2 D(.T-t;t) [4(8- 13t;2 + 5^;4) + (11- 5D2)(a;- 1)^)2] 

^^--^ - 713^ 72 [xi(l - .2) + _ ^1^4d) 

^ VXT' v{l-v') [8x1(1 -v^) + n{x- vtf] 
^ 72[xl{l-v^) + {x-vtYf'^ 

^/A^2 t;(l-t;2)(x-z;t) [5a:i(l-^2)^g^^] 

\/r^ 72[xl{l-v^) + {x-vty^f'^ 

(T ) = -^ v{l-v-){x-vt)xl 

\/r^9[r2(l-i;2) + (x-,;t)2]5/2 

where (x_l,'(?) are polar coordinates for the Xix-g, plane. A strange feature of (164a) 
is that it exhibits a transition from a region of energy depletion behind the quark, 
to a region of energy depletion in front of it as the quarks velocity decreases. When 
> 5/8 there is a buildup of energy density ahead of it, forming a "bulldozer 
effect." See figure|3] As it slows down extra lobe-like features appear until v^ < 5/13 
where the energy buildup is behind the quark, creating an "inverse-bulldozer" effect. 
See figure [4] Recall that the speed of sound in a conformal fluid is = 1/3, so 
that this transition occurs at velocities which are higher than the speed of sound. 
This indicates that the features we are seeing are not hydrodynamic in nature. 
A more detailed analysis of the deviation of the energy density from linearized 
hydrodynamics can be found in the literature.'^ We will see in section [s] that it 
is probably the near field of the stress-energy tensor which dominates high-angle 
emission of hadrons. It would certainly be interesting to understand the physical 
mechanism behind this near-field behavior. 



7. Numerical results for the holographic stress tensor 



Expression (1641 and the Fourier transform of (1241 capture the near-field and far- 
field asymptotics of the stress-energy tensor. In 6.1.2| we have seen an indication 
that far from the moving source the energy-momentum tensor exhibits hydrody- 
namic behavior. In |6.2[ we have seen that the near-field stress tensor exhibits non- 
hydrodynamic behavior with interesting features, like the multi-lobe structure in 
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figure |4] In the intermediate regime, there is a transition region between hydrody- 
namics and whatever short-distance physics governs the near field. To probe this 



region, one needs solutions to (44 1 for values of K where no analytic asymptotic 



treatment is available. We have obtained such solutions numerically. First, (52a I, 



(52b I, (52d), (52el, and (52gl were solved, and Qa, Qd, and Qe were obtained. 



Then, the resulting momentum space stress-energy tensor was passed through an 
FFT to position space using a 128'^ grid. Such a computation has been carried out 
for the energy densitjEMU and for the Poynting vectorHHHl 
Consider the normalized energy density 



£ 



(165) 



where {T^^) is defined in (89 1 as the stress-energy tensor of the system, minus the 
stress-energy of the thermal bath, minus the divergent delta-function contribution 



( 92 1 at the position of the moving quark. It is convenient to decompose this rescaled 



energy density into a Coulombic term, a near-field (large momentum) term, a far- 
field (small momentum) term, and a residual term: 



ilomb 



UV 



(166) 



The Coulombic term represents the contribution coming from the near field of the 
quark: 



Coulomb 



= -2 



247r 



+ (1-^;^)- 



247r 



(167) 



which is what we found in ( 162 1 converted to momentum space. It can be read off 
of the 0{K) terms in (1591. The far-field term, £m, scales like 0{K~^) at large 



momentum and asymptotes to (125a I at small momenta]^ There are many possible 
expressions which satisfy the above criteria. Taking note of the resummation ( 168 1 
we used 



IR 



1 3^7^1(1 



(168) 



where /xir is a typical scale where (125a) stops being valid. We used /iiR = 1. 
Similarly, £uv scales like 0{K^) at small momenta, and asymptotes to (1591 at 



large momenta. To regulate the large momentum expressions in the IR we made 
the replacement 

-n/2 



1 



(if2 _ „2^2)„/2 (i^2 



i/2 



1 



(169) 



where the last term in (169) is expanded to order 0{K °). Similar to (168), /iyv 



is a cutoff scale which we set to 1. After the replacement ( 169 1, the first two terms 



''Actually, we require that it asymptote to | |125a| l up to a momentum-independent constant. 
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in a series expansion of the energy density take the form 



{2 + v^)^K^ 



24 



(170) 

Once fuvi '^'iRi ^^^id fcouiomb are known we can numerically compute E^csi which can 
be fed through a three-dimensional FFT with controllable errors because it is abso- 
lutely integrable. We then add back to the real space numerical expression for f^es 
the real space version of 5ir, 5uv, and f couiomb to obtain the energy density in po- 
sition space. The Fourier transform of the large momentum asymptotic expressions 



can be carried out using (1331. As explained in section 6.1.2 Fourier transforming 



the sound pole is difficult due to the cubic term in the denominator — the term as- 



sociated with the shear viscosity. To convert the the sound pole structure of (168) 
to real space we first rewrote it as 



£ 



IR 



Kl 



m(Ki 



(171) 



then Fourier transformed in the K± direction using ( 133 ), and resorted to numerics 



to FFT the remaining Ki coordinate. This was carried out on a line with 1944 
points with Ki ranging from —20 to 20. In figures [3||5] we show the energy density 
(with the Coulombic field fcouiomb subtracted) at various spatial scales. 
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Fig. 3. (Color online.) A contour plot of the energy density near the moving quark, with the bath 
and the Coulombic contributions subtracted.!^ Hed signifies energies above the background value 
of the plasma while blue signifies energies below the background value of the plasma. The black 
dot at Xi = Xp = marks the location of the quark which is moving at a constant velocity of 
V = 0.58, just above the speed of sound. We work in dimensionless units where Xi = wTxi and 
Xp = TtTx^. The dashed green line shows the presumed location of the Mach cone. 
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E for v=0.75 
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Fig. 4. (Color online.) A contour plot of the energy density near the moving quark, with the bath 
and the Coulombic contributions subtracted i^*? Red signifies energies above the background value 
of the plasma while blue signifies energies below the background value of the plasma. The black 
dot at Xi = Xp = marks the location of the quark which is moving at a constant velocity of 
V = 0.75, well above the speed of sound. We work in dimensionless units where X\ = wTxi and 
Xp = nTx±. The dashed green line shows the presumed location of the Mach cone. 
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Fig. 5. (Color online.) A contour plot of the energy density far from the moving quark, with the 
bath and the Coulombic contributions subtracted.!^ Red signifies energies above the background 
value of the plasma while blue signifies energies below the background value of the plasma. The 
black dot at Xi = Xp = marks the location of the quark which is moving at a constant velocity 
of t! = 0.75, well above the speed of sound. We work in dimensionless units where Xi = nTxi and 
Xp = ttTx^. The dashed green line shows the presumed location of the Mach cone. 
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The components of the energy flux can be treated in a similar manner: we deflne 



and decompose S into 



S — Scou 



loinb 



uv 



S 



IR - 



Sic 



(172) 



(173) 



The Coulombic expression for the Poynting vector is given by the 0{K) terms in 

exprc 



( 160 1 and ( 161 1. The small momentum expressions are given by 



Si 



IR 



Mir 



2v 1 + iKi/Av 2v 



iKi/Av 



(174) 

where we have set /iiR = 1 and by "(regulators)" we mean terms containing the 
regulator /iip, analogous to those in (1681 and (1741. The large momentum ex- 
pressions are given by applying ( |169[ ) to (1601 and (1611. As was the case for the 
energy density, we used /iuv — 1- The real space results for the Poynting vector for 
— 3/4 are shown in figure [6] 
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Fig. 6. (Color online.) Contour plot of the magnitude of the Poynting vector, with the Coulombic 
contribution subtracted.lSH The magnitude of the Poynting vector goes from red (large) to white 
(zero) while the arrows show its direction. The dashed green line shows the presumed location of 
the Mach angle, and the blue line shows the location of the laminar wake — as dictated by its large 
distance asymptotics. 
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8. Hadronization, jet-broadening, and jet-splitting 

There is a significant gap between the results of sections [5}|7] and experimental 
data. Before reviewing recent attempts to bridge this gap, let's briefly summarize 
some of the relevant data. There are of course more authoritative summaries in the 
experimental hteraturelUMIlMIIllMl 

The data seem to reveal a phenomenon of "jet-splitting," whereby an energetic 
parton traversing the medium deposits so much of its energy through high-angle 
emission that — with appropriate momentum cuts and subtractions — the extra par- 
ticle production due to the parton is at a minimum in the direction of its motion, 
and reaches a maximum at an angle roughly 1.2 radians away. Jet-splitting is most 
simply illustrated through two-point histograms of the azimuthal angle separat- 
ing a pair of energetic hadrons close to mid-rapidity. To understand the phenomena 
better, it is useful to examine two landmark studies of these histograms: one from 
STARPand one from PHENIX.E^ 

In the STAR analysis, fairly inclusive momentum cuts were considered: under 
one set of cuts, the less energetic of the two hadrons was required to have trans- 
verse momentum greater than 150 MeV/c. The resulting data show a peak for nearly 
coUinear hadrons that is approximately the same shape for central gold-gold colli- 
sions as for proton-proton collisions: see figure [7j This "near-side jet" feature can 
reasonably be supposed to arise from vacuum fragmentation effects. The two-point 
histogram also shows an "away-side jet" feature around A(/) = tt which is substan- 
tially broader for central gold-gold collisions than for proton-proton. Neither the 
"near-side jet" nor the "away-side jet" are reconstructed jets in the usual sense; 
instead, they are ideas that help explain the main features of histograms assembled 
from millions of events. It is usually assumed that the typical event contributing 
to the histograms involves a hard scattering event where one parton escaped the 
medium without much interaction, producing the highly energetic "trigger hadron," 
and the other parton interacted substantially with the medium before generating 
an "associated hadron" in the vicinity of — tt. The upshot is that with inclu- 
sive momentum cuts, there is substantial broadening of the away-side jet, but not 
jet-splitting: associated hadron production is still maximized, or statistically indis- 
tinguishable from its maximum, at Ai^ = tt. With tighter momentum cuts on the 
associated hadron, the data used in the particular STAR analysis under discussion 
show striking jet-broadening, but the scatter in the data is sufficient to prevent firm 
conclusions from being drawn — from this particular study — about whether there is 
jet-splitting. (Subsequent STAR analyses of both two- and three-point hadron cor- 
relators provide strong evidence in favor of jet-splitting.'^Zl^^Ml) 

The PHENIX analysiplis similar to the STAR analysis,'s31but with more restric- 
tive cuts: in particular, the less energetic hadron was required to have transverse 
momentum greater than 2.5GeV/c. The resulting histograms, with a zero-yield- 
at-minimum (ZYAM) subtraction, show a distinct minimum in associated hadron 
production at At/) = tt, with a broad maximum in the ballpark of A(/) = tt — 1.2. 



Energy loss and the gauge-string duality 45 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 
A (|) (rad) A ()) (rad) 



Fig. 7. (Color online.) Top: The STAR analysiJ^ shows substantial broadening of the away-side 
jet. Reprinted figure 1 with permission from J. Adams et al., Phys. Rev. Lett. 95, 152301 (2005), 
|http://hnk.aps.org/abstract/PRL/v95/pl5230l] Copyright 2005 by the American Physical Soci- 
ety. 

Bottom: The PHENIX analysis^^ shows jet-sphtting for sufficiently central events. Reprinted 
figure 2 with permission from S. S. Adler et al, Phys. Rev. Lett. 97, 052301 (2006), 
|littp://link. aps.org/abstract/PRL/v97/p052301 Copyright 2006 by the American Physical So- 
ciety. 

In reference to the images in this figure, we note in accord with the publisher's guidelines that 
"Readers may view, browse, and/or download material for temporary copying purposes only, pro- 
vided these uses are for noncommercial personal purposes. Except as provided by law, this material 
may not be further reproduced, distributed, transmitted, modified, adapted, performed, displayed, 
published, or sold in whole or part, without prior written permission from the American Physical 
Society." 
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This jet-splitting persists down to roughly 50% centrality, meaning that it occurs for 
events where the impact parameter is less than about 10 fm. The ZYAM subtraction 
is an important part of the analysis, especially for the less central events. The reason 
a subtraction is needed is that for a non-central collision, there is an angular mod- 
ulation of single-particle yields, approximately proportional to 1 -I- W2 cos 2(f), where 
the zero of (j) coincides with the azimuthal direction of the impact parameter, and 
the elliptic flow coefficient V2 depends on the transverse momentum and species. 
Two-point hadron correlators receive a contribution from single-hadron yields. The 
ZYAM scheme is to subtract a multiple of the appropriate product of single-particle 
yields. The multiple is chosen so that the resulting histogram has one bin with zero 
net events, while all other bins have a positive net number of events. 

A natural hypothesis is that the high-angle emission leading to either jet- 
splitting or jet-broadening can be described in terms of a sonic boom in the 
mediumE^iZU Two related difficulties afflict this idea. First, it's hard to get a sonic 
boom with a big enough amplitude to account for the data®^ with reasonable rates 
of energy loss;^^ it should be noted however that not all investigators agree on this 
point and that there are some phenomenological models based on sonic booms 
that flt the dataElH^ Second, one usually finds a diffusion wake with comparable 
strength to the sonic boomEnHHniZll At least in a static medium, it is hard to get 
jet-splitting in the presence of a significant diffusion wake. 

To compute the relative strength of the diffusion wake and sonic boom for the 



heavy quark in the SYM theory, we go back to the conservation equation ( 99 1 , 



F" = izl Hm K„(r^" 



K 



(176) 



In section [6. 1.2| we saw that at small K, the components of the stress-energy tensor 
may be decomposed into terms containing sound poles at ~ SKiv"^ and terms 
associated with a wake which have a pole at Ki ^ 0. If {T^) is the leading, small 
K contribution to the stress-energy tensor, then we may decompose 
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Since only the terms in (1771 contribute to the total drag force in (1761, this gives 



a natural division of the total drag force: 

Konnd = *4 lini /^^(T^rnd) ^wakc = *4 1™ K^{TZlo) , (180) 



From (1241, (125l, (1261, and (1271, we find that 



In our conventions, the zero component of gives us the total rate of change in 
the energy density plus any energy flux going out of the system. Thus, the ratio of 
energy going into sound waves to energy going into the wake is 

l + v^:-l. (182) 

While sound modes carry energy away from the moving quark, the wake feeds 
energy in toward the quark. While this may seem counter-intuitive, in some sense 
it's obvious: the diffusion wake consists of a flow of the medium forward toward 
the quark. The forward-moving momentum in the diffusion wake is the momentum 



deposited by the quark at earlier times. Qualitatively, (182) says that the diffusion 



wake and the sonic boom have comparable strength. When comparing ( 182 1 to the 



scenarios of energy loss in the literaturejUEH one flnds that (182) quantitatively 



matches a scenario where the relative strength of the wake is so large that it washes 
out features of jet-splitting associated with the sonic boom. It may be signiflcant, 
however, that the medium is inflnite and static, both in our work and in the linear 
hydrodynamic scenaricJ^ that our results match onto at large length scales. 

In light of the difficulties in explaining the data with a "boom and wake" model, 
focused on the hydrodynamic regime, it is natural to investigate the effect on hadron 
production of the region of the medium close to the moving quark where hydrody- 
namics is inapplicable. This has been pursued in a series of works 
which we briefly summarize in the next few paragraphs. The first idea is to sub- 
tract away the leading order Coulombic contribution to (Tmn)- Up to an overall 
multiplicative rescaling, these are the quantities we denoted £couiomb and 5couiomb 
in section [7] The justification for this is that these fields describe the energy of the 
energetic parton itself, not the energy lost from it. The remaining energy density, 
which we will denote as Csub, can be split up as Csub = Cbath + Ae. (Note that in 
contrast to our definitions of £ and its variants, Csub explicitly includes the contri- 
bution from the bath.) The Poynting vector Ssuh, with the Coulombic contribution 
subtracted away, is non-zero only because of the presence of the quark. The basic 
plan is to use the Cooper-Frye algorithm*^ to convert string theory predictions for 
(ssub, 'S'sub) into a spectrum of hadrons. 

The Cooper-Frye algorithm is based on converting a fiuid element at tem- 
perature T and with local four-velocity [7™ into hadrons according a Maxwell- 
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Boltzmann distribution in the local rest frame: 

dN 



f{PT,4') 
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oo 



Jq Jo J-oo 



where N is the number of hadrons, and we have set 

P™ = (pt Pt cos(7r — 0) pt sin(7r — 0) O) 

and 



U" 



(U° C/j_cos<y9 C/j_sin(p) 



(183) 



(184) 



(185) 



Note that because of our choice of mostly plus signature, the energy of the hadron in 
the local rest frame of the fluid is —U"^Pm- Also because of this choice of signature, 
we are obliged to include an explicit minus sign in the first integral expression of 



(1831 



To understand (183h(185|, it helps to refer to figure [s] The momentum of 
the associated hadron is P™, and (183) is written in the approximation that the 
associated hadron is massless — an excellent approximation since a typical hadron 
of interest is a pion with ~ 3 GeV / c. The rapidity y is related to the angle from 
the beamline 9 by 



tanh y = cos 9 . 



(186) 



(Note that rapidity y has nothing to do with the depth coordinate y — z/zh used 
in previous sections.) The freeze-out surface is chosen to be a slice of constant a;° 



in (1831. This is the best motivated choice for an infinite, asymptotically static 



medium. In an expanding medium, a more usual choice is a fixed-temperature 
surface with the temperature set close to the QCD scale. For isochronous freeze-out, 
the measure dS'" is simply dx^dx^dx^ (l O), and in passing to the second 



line of (183) we have simply expressed the metric on R'^ in radial coordinates. It 



is important to realize that the azimuthal angle ip around the direction of motion 
of the parton (assumed to be in the -l-a;^ direction, as usual) is different from the 
azimuthal angle <j) around the beamline. We have omitted in ( 183 1 a subtraction of 



the contribution of the bath to hadron production, which depends on pT but not 
on (j). We have also not attempted to normalize f(pT,(t>)- doing so would involve 
partitioning over the spectrum of hadrons. 



The information from string theory enters into (183) in two ways. First, the 
local four-velocity U"^ is the local rest frame of the medium, in which the Poynting 
vector vanishes. Second, the temperature T is the temperature in this local rest 
frame, deduced by plugging the energy density T^^V^U^ into the equation of 
state. Evidently, one needs all components of T™^ in order to precisely determine 
t/™ and T. This is a problem since only the m — Q row of the stress tensor has been 
computed in full.^^ Another problem is that close to the quark, the medium is 
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Fig. 8. (Color online.) A sketch of the coordinate systems used to describe away-side hadron 
production. If the trigger hadron is at = 0, then a helpful relation at mid-rapidity (9 = 7r/2) is 
= ■7T — <p, where A<p is the azimuthal separation between the trigger and associated hadrons. 



presumably far from equilibrium, so using Cooper-Frye seems somewhat perilous. 
We will return to a discussion of these two issues below. 

To understand how hydrodynamical and non-hydrodynamical effects contribute 
to the spectrum of produced hadrons, one must have some notion of where 
the boundary is between hydrodynamical and non-hydrodynamical regimes. This 
boundary is presumably not sharp. Three considerations have gone into identifying 
an appropriate boundary: 

(1) The non-hydrodynamical region can be chosen as the region where Ae/etath is 
less than some constant of order unity. For 

A = 5.5, iV = 3, rsYM = 200MeV, (187) 

a preferred choice is 

Ae/ebath < 0.3 . (188) 



Here and below, we will describe the region defined by ( 188 1 with the parameter 



choices ( 187[ ) as the "Neck." It extends roughly over —1 < Xi < 0.5 and 



< < 1.7. (For TsYM = 200 MeV, X = 1 corresponds to a; = I/ttTsym 
0.31 fm.) 
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(2) The Neck region can be compared with the region where the Knudsen number 
exceeds some constant of order unity. An appropriate version of the Knudsen 
number in the current context is 

Kn = r '^:'^^"'^' ■ (189) 

I 'S'sub I 

Here the sound attenuation length F is the same as the one discussed following 



(120): r = Arj/SsT — l/SnT. For the choice of parameters (1871, examination 



of the near- field expressions ( |164 l shows that the region where Kn > 1/3 is 



somewhat bigger in the Xi direction than the Neck. However, corrections to 
the near- field approximation to S'sub may not be negligible for Xi and/or X± 
of order unity. 

(3) The Neck region can be compared with the region where the constitutive rela- 
tions of hydrodynamics break down. Given all components of the stress 
tensor, with the Coulomb field subtracted away, there is a straightforward pro- 
cedure for testing the constitutive relations. First determine the local velocity 
field t/™ by passing to the local rest frame of the fluid. Let (T'"")^ be the 
subtracted stress-energy tensor in the local rest frame. The energy density is 
read off immediately as (T^")^; the pressure is deduced from the equation of 
state; and the shear viscosity contribution to the space-space parts of the stress 
tensor can be determined from [/™ and its gradient. Deviations of the space- 
space components of (T"*")/^ from the combined contribution of pressure and 
shear viscosity are measures of the failure of hydrodynamics. 
A studjl^ of the near-field expressions (|164[) for a somewhat different choice 



of parameters from (1871 (namely v — 0.99, A = 'Sn, N — 3) concludes that 
deviations from hydrodynamics are appreciable out as far as X ^ 8. However, 
the near-field expressions definitely cannot be trusted at such large distances^ 
This analysis could therefore be considerably improved if all components of 
-^sub were computed directly from string theory. 

The main conclusion to draw from points 2 and 3 is that in the Neck region, 
the subtracted stress tensor is essentially unrelated to hydrodynamics. Instead, the 
physics may be presumed to be dominated by strong coherent color fields combined 
with responses of the medium to strong field gradients: hence the term "chromo- 
viscous neck." 

To return to hadronization: The Cooper-Frye integral over R'^ can be split 
into the Neck region and the "Mach" region — which is everything else. In the Mach 
region, where the energy density comes mainly from the bath, a good approximation 



'A computation of the energy radiated from a moving quarlfiS shows that agreement with hn- 
earized hydrodynamics is already fairly good at X ~ 5. Based on results of an earlier study,[^the 
onset of reasonable agreement with linearized hydro occurs near X = 3. 
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to the local rest frame can be found by setting 

U=^^. (190) 

4 Ebath 

In the neck region, this approximation is less reliable, but because space-space 
components of J"™" are not available from a string theory calculation, it is hard to 
give a better motivated prescription for determining the local rest frame. With the 



choice (1901, the result is that the Neck contribution to the Cooper- Frye integral 
leads to a distinctive double-peaked structure in /(pt, for pT- This is remarkable 
when compared to the single-peaked structure emerging from a computation in a 
perturbative QCD framework based on Joule heating,!^ which is similarly passed 
through the Cooper- Frye hadronization algorithm. See figure |9] 

The double-peaked structure from the Neck region of the trailing string stress 
tensor has nothing to do with the Mach cone. It doesn't occur at the same angle: 
for example, at w = 0.58, the Mach angle 1) = cos~^ Cg/v is very nearly zero, but the 
Neck region still produces a double peak structure (not shown in figure |9| about a 
radian away from (j) = tt. When v is very close to 1, the double peaks get closer to 
(j) = TV. This is reminiscent of the structure observed at large K in Fourier spacejUEH 
but the peaks observed in the predicted hadron spectra are more widely separated 
than the ones in Fourier space. 

While the hadronization studieij^MZr^r^lsoisilEMI give valuable insight into the 
relation of the trailing string to high-angle hadron emission from an energetic par- 
ton, it is not claimed that the results are fully realistic, or that a direct compar- 
ison to untagged dihadron histograms, like the ones in figure [7} is justified. Let 
us review the potential difficulties. First, the trailing string describes an infinitely 
massive quark that propagates at a constant velocity through an infinite, static, 
thermal medium. Fluctuations leading to stochastic motions of finite mass quarks 
may significantly affect the results. Also, there may be an effective "form factor" 
for massive quarks that partially smears out the field close to the quark. In addi- 
tion, it is not obvious that Cooper-Prye is justified, because the crucial effect comes 
from the non-equilibrium part of the medium (the Neck). Subtracting away the 
Coulomb field is certainly well-motivated physically, but it is possible to maintain 
some skepticism about whether it is the correct prescription in combination with 



Cooper-Frye. Finally, the approximation (1901 to the local rest frame is imprecise 
in the Neck region. 

Despite these potentially serious issues, the punchline of the phenomenological 
studiejl ^ l ^^ l ^^ l ^^ l ^° l ^^ l ^^ l ^'^ ] seems to us likely to be robust: the near-quark region has 
a substantially greater tendency toward high angle emission in the trailing string 
treatment than in the perturbative QCD treatment based on Joule heating. Modulo 
concerns already expressed, the near-field contribution to high-angle emission is 
stronger than the contribution of the hydrodynamical regime, and it results in a 
significant double-hump structure, reminiscent of jet-splitting. 
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Fig. 9. (Color online.) Hadron production based on a Cooper-Frye hadronization of a perturba- 
tive QCD calculation and of trailing string results (AdS/CFT).EIl Bath contributions have been 
subtracted away, and the curves have all been normalized to have the same maximum when the 
Cooper-Frye integral is carried out over the entire volume accessible to each computation. The 
Neck region for the perturbative QCD calculation is again defined as a region close to the quark 
where deviations from hydrodynamics are significant. 



9. Conclusions 

Let us conclude by addressing the four main questions we raised in the introduction: 

(1) What is the rate of energy loss from an energetic probe? 

For a heavy quark moving at a velocity v, the drag force is -Fdrag = 
_7ry52i2_ii_ fpjjjg jg explained in more detail in section jZj 

(2) What is the hydrodynamical response far from the energetic probe? 

There is a sonic boom with Mach angle i9 = cos~^ Cs/v, where Cs = is the 

speed of sound dictated by conformal invariance. There is also a diffusion wake 
of comparable strength. These points are explained in section [6. 1 .2| and [T] 
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(3) What gauge-invariant information can be extracted using the gauge-string du- 
ality about the non-hydrodynamic region near the probe? 

The expectation values (T"'''^) of the energy density and the Poynting vector 
of the gauge theory stress tensor have been computed with uniformly good 
accuracy across all length scales for several values of the velocity of the heavy 
quark, as we review in section |7] Analytic approximations to the space-space 



components of (T™") are also available at small length scales: see section 6.2 



(4) Do the rate and pattern of energy loss have some meaningful connection to 
heavy ion phenomenology? 

A suitable translation of parameters from SYM to QCD results in estimates of 
energy loss for c and b quarks which are not far from realistic, or which may 
be fully realistic. We summarize these estimates in section |4] 
Studies of hadronization starting from the string theory predictions for the 
energy density and Poynting vector indicate that the trailing string leads to 
significant high-angle emission from the Neck region, close to the quark, sug- 
gestive of jet-sphtting. We describe these studies in section [8] Although it is 
premature to make detailed comparisons to data, it is clearly worthwhile to 
extend and refine both the string theory analysis and the phenomenological 
studies. 
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Appendix A. Notation 

In this appendix we present short explanations of some of the nomenclature and 
mathematical notations used in the main text. 

x^^: The five spacetime coordinates of AdSs-Schwarzschild, usually {t, .x^ , z). 
Pm'- The four-vectors for position and momentum in R'^'^. We use mostly plus sig- 
nature, so (for example) vi"^"'PmPn = —E"^ + ]f = — m^. 
X, p: The three-vectors for position and momentum: spatial components of x™ and 

Xi: This could mean either Gif^x^^ — ^x^ or rji^x'^ — x^. Our convention is to 
prefer the latter; likewise X2 = x^ and X3 — x^. 

x±: The radial distance from the quark in the X2, X3 plane: xj_ ~ ^ x\~\- x\. Occa- 
sionally we consider the two-vector x\^ = (x2-,x-i). 
z: This is the depth coordinate in AdSs or AdSs-Schwarzschild which is at the 
boundary and has dimensions of length. The AdSs metric is ds^ = ^{—dt^ + 
dx^ + dz'^). 

zh'- The depth of the horizon in AdSs-Schwarzschild, related to the temperature by 
T = I/ttzh- 

y: Usually, a rescaled depth coordinate in AdSs-Schwarzschild, defined hy y = 
z/zh- But in section [s] we use y to denote rapidity, i.e. tanhy = Pz/E where 
Pz is the momentum along the beampipe and E is the energy, 
r: We use r to indicate a radial separation in R^. Some authors use r to denote 
the depth coordinate r = L^/z in AdSs. 
AdSs: Five-dimensional anti-de Sitter space, the maximally symmetric negatively 

curved spacetime in 4 -I- 1 dimensions. Its metric is given by ([s]) with h — 1. 
SYM: An abbreviation for "A/" = 4 super- Yang-Mills theory in four dimensions," which 
is the theory controlling the low-energy excitations of D3-branes. 
TV: Usually, the number of colors: = 3 in QCD. An exception is that in section 8, 
we use N to indicate the number of hadrons predicted by the Cooper-Frye 
algorithm. 

gvM- The gauge couphng of SYM, normalized so that g^y^N — /a'"^, where N is 
the number of colors. 

The gauge couphng of QCD. We also use Ug = gl/^T^- 
The 't Hooft couphng, A = ^ym-^- 
The radius of curvature of AdSs. 



A 
L 
a' 

K 

R 



The Regge slope parameter of fundamental strings, see (14 1. 

The five-dimensional gravitational coupling. 

The spacetime metric of AdSs or AdSs-Schwarzschild. 

The Ricci scalar in AdSs or AdSs-Schwarzschild. We also use the Ricci tensor 
R^i, and the Riemann tensor RafifSu- Our conventions are R ~ G'^'^R^i, and 
Rfif = G°'^ Rafifju, with signs arranged so that R^i, — —-^G^u in AdSs of 
radius L. 

The "blackening function" for AdSs-Schwarzschild, whose metric is given in 
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T: 

9a(3 
V 

a;* (cr): 



R 



AA- 



'Viz- 
Axial gauge: 



A: 



ay. 
K: 



(|3|. It is given by h{z) = 1 — z'^/zj^. We sometimes think of /i as a function 
of the depth z, and sometimes as a function of y z/zh- h! always means 
h'{y) = ~Ay\ 

The temperature in the dual field theory, which is the same as the Hawking 
temperature of the dual black hole background. The temperature of the AdSs- 
Schwarzschild background ([s]) is T = I/ttzr- 
The worldsheet metric of a string. 
Coordinates on the string worldsheet. 

The speed of a moving quark. Usually we take this motion to be in the -\-xi 
direction. 

Gives the shape of the string that describes the quark in the five-dimensional 
geometry. See (13 1 and (20 1. 

The embedding function for a classical string in AdSs-Schwarzschild. When no 
ambiguity is possible, we denote this embedding function more simply as a;^(cr). 
The nuclear modification factor, defined as the number of particles produced 
(usually at a particular value of pt and in a specified range of rapidity) in a 
collision of two nuclei with atomic number A, divided by the number produced 
in a proton-proton collision scaled up by the effective number of binary nucleon- 
nucleon collisions in the heavy-ion collision. 

The five-dimensional stress-energy tensor, not to be confused with the expecta- 
tion value of the boundary stress-energy tensor (Tmn)- For the trailing string, 
T^^ is given in ([42 1. 

The Fourier components of r^jy. See (46 1 and (48 1. 

Small metric perturbations around AdSs-Schwarzschild. 

The Fourier components of h^^^, defined by analogy with (46 1. 

A gauge choice for the metric perturbations where /i^^ = 0. 

The Fourier components of h^^^ in axial gauge, up to a normalization factor. 

See (49). We think of the Hmn as functions oi y = z/ zh , not of z. 

The even tensor mode combination of metric perturbations in axial gauge: see 

(51a). Similar definitions for Bi, C, Di, and Ei follow, and gauge-invariant 



combinations B, D, and E can be found in (61l-(62l. These are all functions 

of J/ = z/zH- 

The scalar master field. A master field is a gauge-invariant combination of 
metric fiuctuations in AdSs-Schwarzschild with a simple equation of motion. 
We encountered four other master fields: V'y'™, V'y^'^i ''Pj 

1/ (27ra'\/r 



and tf}!^'^'^: see ((67 1. 



A recurring normalization factor given by a^, = i/ {^'Zna' \l I — 
The momentum conjugate to x/ z^; used as a dimensionless wave- number to 
parameterize the three-dimensional Fourier space used to describe the medium's 
response to the quark. 

The magnitude of the component of K perpendicular to the motion of the 
quark. 

The asymptotic values of Hmn aX y — Q. The boundary condition Rmn — says 
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Prr 



( T^hydro \ 

\-'- jmn- 

\ mn/ 
/n ■ 

Rx- 

Px- 
Qx- 

Pt- 
ip: 



I,K: 

Neck: 



that the four-dimensional metric which the boundary gauge theory experiences 
is flat Minkowski space. 

The coefficients of in a small y expansion of Hmn ■ See ( 76 ) and ( 77 1 . They 



are related to the divergent contribution to the boundary stress-energy tensor 
given in ( 90 1 . 

Qmn- The coefficients of y^ in a small y expansion of Hmn when Rmn ~ 0. They are 



related to the expectation value of the boundary stress-energy tensor by (91 1 



The metric of the boundary conformal field theory, usually set equal to the 
Minkowski metric rfmn with mostly plus signature. 

The four-velocity of a heavy quark moving through the thermal medium. 
The one-point function of the stress-energy tensor in SYM in the presence of 
the moving quark. We find it convenient to decompose it into three pieces given 



in (891 



A stress tensor which satisfies the hydrodynamic constitutive relations, (1201. 
^jihydro.SYM-j^^^ jg ^-j-^g hydrodyuamic contribution to the stress tensor of the 
SYM theory. 

The Fourier modes of the contribution of the moving quark to the stress-energy 
tensor in the dual field theory minus the divergent piece corresponding to the 
infinite mass of the quark. See (89). It can be computed from Qmn through 



(95 



The source term for the energy- momentum tensor, iK^T™^ — /". Various 
superscripts specify which contribution of the energy-momentum tensor is being 
sourced. For example, f^y^"^" sources (T^yd™)^^^ 

The coefficient of the leading homogeneous solution for various linear combina- 
tions of Hmn and their derivatives, such as A, B, Bi, ipx, etc. Our boundary 
condition is Rx — 0. 

The analog of Pmn for various linear combinations of Hmn and their derivatives. 
The analog of Qmn for various linear combinations of Hmn and their derivatives. 



See (1021 and (108). Also, we use the notations Qd — Qdi and Qe ~ Qei- 
The component of momentum perpendicular to the beamline. 
An angular coordinate in momentum space, sinz9 = Ki/K, or in real space, 
sin — xi/x. 

The azimuthal angle around the direction of motion of an energetic quark. Our 
usual convention is that the energetic quark moves in the -l-a;^ direction, and 
then tan ip — / K2 or X3 /x2 . 

The angle of a trajectory relative to the beam. = 7r/2 is mid-rapidity. 

The azimuthal angle around the beam. The angular variable in dihadron 

histograms is the separation in </) between two hadrons. 

Modified Bessel functions of the first and second kind. J is a Bessel function of 
the first kind, and L is a modified Struve function. 

The neck is the region near a moving quark where the response of the medium is 



non-hydrodynamical. In practice, for the choice of parameters (187), the Neck 



can be defined, as in (188), as the region where the energy density, excluding 
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the Coulombic contribution, exceeds 1.3 times the asymptotic energy density 
of an infinite static bath. 
[7™: The four- velocity of a fluid element, usually defined so that it vanishes in pre- 
cisely the same Lorentz frame in which the Poynting vector vanishes. 
S: The energy density and Poynting vector, rescaled to make them dimensionless, 
with contributions from the thermal bath excluded: see (165) and (1721. £ can 
be decomposed into a sum of contributions from the Coulomb field of the quark, 
subleading UV effects, IR effects, and a residual quantity fjcs, as in (1661. An 
analogous decomposition can be performed on S. 
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